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FOREWORD 


This report was prepared by Wyle Laboratories/ Research Division/ Huntsville, Alabama for 
the Unsteady Gasdynamics Branch, Aero-Astrodynamics Laboratory, Marshall Space Flight 
Center (MSFC), National Aeronautics and Space Administration (NASA), Huntsville, Alabama 
under Contract NAS8-21026. The work was performed under the direction of Messrs. G. 

Wilhold - Deputy Chief of the MSFC Unsteady Gasdynamics Branch, and J . Jones - Technical 
Monitor of the contract. 

This report contains the results of a wind tunnel test conducted at the Arnold Engineering 
Development Center (AEDC) during the period from May 22 through J une 7, 1_?68 . The tests 
were conducted under the direction of Mr. H. C. Dubose of ARO, Inc., contract operator of 
At DC, with Mr. T. R. Brice as the project engineer. 

The author wishes to express his appreciation and to extend thanks to all persons at NASA-MSFC, 
AEDC, and Wyle Laboratories who participated in this program of study. Those particularly 
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J. Black of ARO, Inc. - AEDC, and Messrs. A. Jolly, D. Wallis, J. Matzkiw and S. Dendrinos 
of Wyle Laboratories. 
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ABSTRACT 


A wind tunnel investigation was conducted to study the flow field induced by three-dimensional 
protuberances at transonic Mach numbers. The perturbed flow environments of both the pro- 
tuberances and the surrounding structure were analyzed utilizing static- and fluctuating-pressure 
measurements and oil-flow visualization techniques. Test results for the generalized configura- 
tions which consisted of 2- , 4-, and 8-inch diameter cylindrical protuberances are discussed 
herein. The range of test variables consisted of protuberance heights from zero to 2 diameters, 
free-stream Mach numbers from 0.60 to 1 .60, and unit Reynolds numbers from 1 .5 to 4.5 million 
per foot. The static pressure results are presented in coefficient form to show the axial pressure 
distributions associated with the protuberance induced flow field as well as the distribution of 
surface pressures over the wall of the protuberances. Fluctuating pressure measurements are 
presented in the form of overall RMS fluctuating pressure coefficients, power spectra, cross- 
power spectra, and narrow band convection velocities. From the static- and fluctuating- 
pressure measurements and the oil-flow patterns, the structure of the perturbed flow field has 
been defined. One of the most significant findings is that the upstream separated flow field 
induced by three-dimensional protuberances consists of a complex, multiple vortex system which 
generates fluctuating pressures an order of magnitude greater than those observed in two- 
dimensional separated flows. The extreme fluctuating pressures encountered within the three- 
dimensional separated flow field are attributed to the shear interaction of two major vortices 
within the separated region. 
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INTRODUCTION 


1 Background 

Protuberances have an infinite number of possible geometries so that a general 
discussion of the flow field is difficult. For example, the Saturn V alone has 
approximately 75 individual protuberances which include reaction control rockets, 
auxiliary propulsion systems, vents, tunnels, etc. (Reference 1). Each protuberance 
generates its own flow field which may interact with the external flow field already 
present and lead to the imposition of large steady and fluctuating loads on both the 
protuberance and the surrounding structure. It is obviously impractical to investigate 
the environment of every protuberance which may be attached to the external surface 
of a launch vehicle, especially if some general features of the flow can be defined 
from a systematic study involving generalized protuberance geometries. In view of 
the free interaction hypothesis first advanced by Chapman, Kuehn, and Larson (Refer- 
ence 2), it appears that the use of generalized protuberance shapes is a reasonable 
approach in studying protuberance flow fields. Essentially, they suggest that boundary 
layers undergo separation in a manner which is independent of the original cause. 

The separation commences well ahead of the protuberance, so that it is reasonable to 
suppose that the flow near the separation point is unaffected by the detail geometry 
of the protuberance. This phenomenon has been proven in experiments for two- 
dimensional supersonic flow and it is expected to apply to the three-dimensional case 
as well . However, even with this simplification, the flow close to the protuberance 
must be expected to be a function of the detailed protuberance shape. It follows that 
a definitive study of protuberance flow would involve measurements on the protuberance 
as well as on the structure which is in close proximity to the protuberance. 

For the past two years, Wyle Laboratories have been engaged in a research program 
under Contract NAS8-21026 with NASA-MSFC to investigate the steady and unsteady 
aerodynamic characteristics of the flow around three-dimensional protuberances which 
extend into, and well beyond, attached turbulent boundary layers. The first task of 
this study was to conduct an extensive search of available literature for the purpose of 
accumulating data for a variety of protuberance geometries. Results of the literature 
survey indicated that considerable data were available for evaluating the mean flow 
around protuberances at supersonic speeds (see References 3-16, for examples). 

However, little data were available on the fluctuating air loads which may be 
experienced in the vicinity of the protuberances. Further, there appeared to be a 
complete lack of both steady and unsteady aerodynamic data for protuberance flows at 
transonic speeds, particularly for three-dimensional protuberances. This lack of tran- 
sonic data poses a serious problem in predicting the steady and fluctuating air loads in 
the vicinity of protuberances for a launch vehicle, since, in the transonic flow 
regime, a launch vehicle experiences maximum dynamic pressure and, consequently, 
maximum aerodynamic loads. However, for the literature survey, some general 
features of protuberance flows were discovered and the results of this study are 
summarized in Reference 17. 
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Following the survey and analysis of the data from published literature, a small scale 
wind tunnel test program was planned and subsequently conducted in the MSFC 14 by 
14 inch Transonic Wind Tunnel . This test was conducted with three primary objec- 
tives. The first objective was to investigate general features of protuberance flows at 
transonic speeds using generalized protuberance geometries. Right circular cylinders 
of various heights and diameters which project into the flow with the cylinder axis 
normal to the flow direction were selected for this study. The second objective of the 
small scale test was to investigate similarities between the flow fields for generalized 
protuberances and the flow fields for protuberances with specific geometries. For 
this phase of the study, models of the reaction control system (RCS) and the auxiliary 
propulsion system (APS) on the Saturn V were tested at the same conditions as for the 
generalized protuberance study. The third objective was to obtain base data that 
could be used to plan a large scale investigation of protuberance flows. The test 
results for the small scale wind tunnel study are presented in Reference 18. At best, 
the small scale test results revealed only cursory features of the protuberance flow 
fields due to the limitations imposed by the model size on instrumentation density and 
measurement accuracy. This is particularly true for fluctuating pressure measurements 
because of the relatively large size of the microphones. 

.2 Present Test Program 

The large scale wind tunnel test program, reported herein, was planned to obtain more 
detailed static and fluctuating pressure measurements for protuberance flows. For this 
test, an existing model of the same general design as that tested in the small scale 
study, but of considerably larger scale was used. The test was conducted in the 
Propulsion Wind Tunnel, Transonic (16T) at the Arnold Engineering Development 
Center (AEDC), Arnold Air Force Station, Tennessee, from May 22 through June 7, 
1968. 

The wind tunnel test reported herein was based to some degree on the free interaction 
hypothesis discussed in Section 1.1. A basic detailed study was conducted using 
cylindrical protuberances of various heights and diameters for which some mean flow 
data were available. A basic reference length in the separated boundary layer generated 
by a protuberance is the undisturbed boundary layer thickness. Thus, protuberance 
height and diameter must be scaled with boundary layer thickness, and to accomplish 
this the test configuration included protuberances lying both within and well beyond 
the boundary layer. In addition, protuberances having more specific geometries were 
tested which consisted of scale models of the RCS and APS protuberances on the 
Saturn V. 

The primary objective of the test was to determine the steady and fluctuating pressure 
characteristics of the protuberance flow field so that the environments of both the pro- 
tuberances and the surrounding structure could be adequately defined. A compre- 
hensive description of the test program was presented in a pre-test report. Reference 
19. The test was performed according to the test plan with only minor modifications. 
Pertinent excerpts from this reference are included herein for the purpose of clarity 
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and coherence. The protuberances were mounted to a curved panel model which was 
supported above the floor of the wind tunnel test section in a splitter-plate type 
installation. The surface panel and certain protuberances were instrumented with 
static pressure orifices and microphones to obtain the desired steady and fluctuating 
pressure measurements. Also, pitot pressure rakes were installed in the test panel to 
obtain descriptive boundary layer profile measurements. A comprehensive discussion 
of the test apparatus is presented in Section 2.0. The protuberances were tested over 
a range of Reynolds number and Mach numbers. The detailed description of the test 
conditions and procedure is presented in Section 3.0. The data reduction procedure, 
which consisted of both on-line and off-line phases, and an analysis of the precision 
of measurements are presented in Sections 4.0 and 5.0, respectively. A general dis- 
cussion of the test results is presented in Section 6.0. It should be noted that over 
300,000 static pressure measurements and over 100,000 microphone recordings were 
obtained during this test. Thus, any reduction and analysis of the test results must be 
selective. For the present phase of the study, emphasis was placed on the reduction 
and analysis of data obtained from the longitudinal centerline distribution of instru- 
ments which were located in the plane of symmetry of the perturbed flow field. An 
analysis of these data, which defines the salient characteristics of the static and 
fluctuating pressures associated with the protuberance flow field, is also presented in 
Section 6.0. From these data, the basic features of the flow field induced by three- 
dimensional protuberances and similarities and differences between two- and three- 
dimensional protuberance flows have been defined. Finally, conclusions resulting 
from this study and a description of areas to be studied in the future are presented in 
Section 7.0. 
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2.0 


TEST APPARATUS 


2 . 1 Test Facility 

The test was conducted in the Propulsion Wind Tunnel, Transonic (16T) at AEDC. 

Tunnel 16T is a variable density wind tunnel with a 16 ft square test section. The 
walls of the test section are composed of perforated plates which allow continuous 
operation through the Mach number range from 0.55 to 1 .60 with minimum wall 
interference. Details of the test section showing the model location are presented 
in Figures 1 and 2. A more extensive description of the tunnel is given in the AEDC 
Test Facilities Handbook and the latest calibration results are presented in Refer- 
ences 20 and 21 . 

2 .2 Test Article 

The basic test article consisted of a 30-degree segment of a cylindrical shell having 
a radius of curvature of 130 inches. The curved panel was mounted in a splitter-plate 
type fixture which was supported above the Hmnel floor (see Figure 1). The pro- 
tuberances and instrumentation were mounted on a secondary panel (hereafter referred 
to as the test panel) which was 36.00 inches wide and 75.18 inches long and flush 
fitted into the basic panel structure (Figure 3) . The right hand side of the test panel 
(when viewed looking upstream) was attached to the basic structure through hinges 
which allowed the test panel to be pivoted open and, thus, provide easy access to the 
instrumentation and protuberance drive system without having to completely disassemble 
the test article. Also the instrumentation and drive system were attached to the test 
panel so that they were exposed when the model was open as shown in Figure 4. This 
feature of the model design greatly simplified model configuration changes and the 
replacement of instrumentation. Protuberances for the basic study were right cylinders 
of 2-, 4- and 8-inch diameter and were telescoped into the flow over the test panel 
with the axes of the cylinders normal to the flow direction. Geometric details of the 
protuberances are shown in Figure 5. All the protuberances were raised and lowered 
with the same drive system. The drive system, shown in Figures 4 and 6, enabled the 
protuberances to be remotely raised and lowered to any height ranging up to 8 inches. 
Space limitations inside the splitter-plate limited the maximum travel of the drive 
system to 8 inches and since it was desirable to test all of the protuberances at heights 
ranging up to at least two diameters, it was necessary to add an 8-inch extension to 
the basic 8-inch diameter protuberance for certain runs. In the fully retracted posi- 
tion, the top surface of all but the 8-inch diameter protuberance with the 8-inch 
extension were flush with the test panel surface. 

The drive system (Figure 6) consisted of three jack screws and a lift plate to which 
the various cylindrical protuberances were mounted. The drive system was supported 
between the underside of the test panel and a base plate by six support rods and two 
spacer bars. The jack screws were driven by a common 1/3 horsepower, 30 rpm 
electric motor through a chain and sprocket drive system. An O-ring seal was inserted 
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in the test pane! around the protuberance to seal the splitter-plate cavity from the 
top of the test panel and thus prevent leakage of air between the protuberance and 
test panel during the test. 

Protuberance configuration changes were accomplished by removing the protuberance 
cap (see Figure 6) and detaching the base of the protuberance from the lift plate. 

Hole patterns in the lift plate provided for mounting all of the cylindrical protuberances 
on a common centerline as well as rotating the 8-inch diameter protuberance from 
0 degree to polar angles of 15 and 60 degrees (counter clockwise when viewed from 
the top). For the 2- and 4-inch diameter protuberances, the gap between the test 
panel and the protuberance was filled with an instrumented insert panel . The insert 
panels were attached to the test panel in the place of the removable ring (see Figure 6). 
O-ring seals similar to that used for the 8-inch diameter protuberance, were inserted 
between the insert panels and the smaller diameter protuberances to sea! the gap. 

In addition to the basic study, models of the RCS and APS protuberances on the 
Saturn V of approximately 20 percent scale were tested. These models were mounted 
to the same test panel used for the basic study, so that existing instrumentation 
could be utilized. Details of the RCS and APS protuberances are not presented herein, 
since test results for these configurations will be forthcoming in a later report. 

2.3 Instrumentation for Data Acquisition 

Instrumentation in a given configuration of the test article consisted of a maximum 
of 296 static pressure orifices and 122 piezo-electric microphones which were flush 
mounted in the test panel and the 8-inch diameter protuberance. The basic 8-inch 
diameter protuberance and the 8-inch diameter protuberance with the extension are 
the only protuberances that contained instrumentation. Details of the instrumentation 
in the test panel, the insert panels and the 8-inch diameter protuberance configurations 
are shown in Figures 7, 8 and 9, respectively. The test panel contained 216 static 
pressure orifices that were distributed in a polar array centered at the protuberance 
center (Model Station =0), Figure 7a. The polar angles for the static pressure orifices 
are denoted by 0 + and are measured in the clockwise direction when viewing the test 
panel from above. Most of the static pressure orifices were distributed along the 
centerline of the test panel (0 + = 0, and 180 degrees) and to the right side when 
viewed looking upstream. One ray of static pressure orifices was located at 0 + = 270 
degrees so that symmetry of the flow could be determined. The locations shown in 
Figure 7a are nominal dimensions and exact locations are presented in Reference 19. 
Also, the test panel contained 107 microphones that were distributed in a polar array 
along the panel centerline and to the left side of the test panel when viewed looking 
upstream (Figure 7b) . The polar angles for the microphones are denoted by 0” and 
are measured in the counterclockwise direction when viewing the test panel from 
above. This array of steady and fluctuating pressure instruments provided optimum 
utilization of the instruments since the protuberance flow fields were symmetrical 
about the panel centerline. 
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The insert panels that were used for the 2- and 4-inch diameter protuberances 
contained static pressure orifices and microphones that were distributed in polar 
arrays similar to those of the test panel (see Figure 8) . Each insert panel contained 
39 static pressure orifices and 15 microphones. Since the 2- and 4-inch diameter 
protuberances do not contain instrumentation, there were a total of 255 active 
static pressure orifices and 122 active microphones for these two configurations. The 
total number of microphone channels available with the recorder system was limited 
to 122. Thus, microphone distributions in the insert panels provided only qualitative 
information of the unsteady pressure flow field in close proximity to 2- and 4-inch 
diameter cylindrical protuberances. However, detailed definition of the fluctuating 
pressure flow field near the cylindrical protuberances was obtained for the 8-inch 
diameter protuberance which had a relative high density of microphones in close 
proximity to the protuberance. 

Both static and fluctuating pressure instruments were located in the basic 8-inch 
diameter protuberance and the extension of this protuberance. Details of the instru- 
ment locations for the 8-inch diameter protuberance configurations are shown in 
Figure 9. The basic 8-inch protuberance contained 24 static pressure orifices and 6 
microphones distributed at polar angles of 0, 30 and 120 degrees in the wall of the 
protuberance and 32 static pressure orifices and 9 microphones distributed in a polar 
array in the protuberance cap. Thus, for this configuration, there were a total of 
272 active static pressure orifices and 122 active microphones. For the 8-inch 
diameter protuberance configuration with the extension there were 48 static pressure 
orifices and 15 microphones distributed at polar angles of 0, 30 and 120 degrees in 
the wall of the protuberance. All of the microphones for this configuration were 
located in the walls of the extension portion of the protuberance with the base of the 
protuberance and the protuberance cap containing only static pressure instrumentation. 
The protuberance cap contained 32 static pressure orifices as for the case of the basic 
8-inch diameter protuberance configuration. Thus, for the 8-inch diameter protuber- 
ance configuration with the extension there were a total of 296 active static pressure 
orifices and 122 active microphones. Both 8-inch diameter protuberance configurations 
were rotated about the centerline from 0 degree to polar angles of 15 and 60 degrees 
such that static and fluctuating pressures were recorded on the wall of the protuberances 
at polar angles of 0, 15, 30, 45, 60, 90, 120, 165 and 180 degrees. Because of the 
slight curvature of the protuberance cap, it was necessary to maintain the cap at the 
0-degree position. That is, the cap was not rotated with the protuberance. 

Other instrumentation on the test panel consisted of two pitot pressure rakes which 
were extended through the boundary layer for measurement of the boundary layer 
velocity profiles. Details of the rakes are shown in Figure 10. The forward rake was 
remotely controlled so that it could be extended and retracted as desired during the 
test. This rake was designed such that, when in the fully retracted position, the top 
surface of the test panel remained smooth. Further, this rake contained 15 pitot pres- 
sure probes distributed above the centerline of the test panel as denoted in Figure 10. 
The pressure for the static orifice in the test panel immediately aft of the forward 
rake was recorded with the rake in the retracted position, and was used together with 
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the pitot pressures in the computation of the boundary layer velocity profiles in the 
upstream vicinity of the test panel . The aft rake was a fixed device and was present 
throughout the test. The aft rake contained 25 pitot pressure probes which were dis- 
tributed above the panel centerline as shown in Figure 10. The pressure for the most 
aft static orifice in the test panel was used in conjunction with the pitot pressures to 
compute the boundary layer velocity profiles for the aft region of the test panel. 

A pitot-static probe was extended into the secondary flow between the splitter-plate 
and tunnel floor to monitor the flow in this region. The previous small scale test 
indicated that the secondary flow may be partially blocked at Mach numbers near 
1.0. By monitoring the pitot-static probe beneath the splitter-plate and the boundary 
layer rakes above the test panel, the effect of any flow spillage over the leading edge 
of the model could be readily determined. The static and total pressures (forward 
rake only) were connected to a system of pressure scanning switches (Scani-valves) 
which were mounted to the underside of the test panel . A total of ten 48-port Scani- 
valves were available to record the pressures. A typical static pressure installation is 
shown in Figure 11a. The total pressures for the aft rake were conditioned through the 
AEDC-PWT Precision Pressure Balance (PPB) system. This system, which was considered 
to be more accurate than the Scani-valve system, was also used to compute a reference 
pressure which was compared to its counterpart on each of the Scani-valves. The out- 
puts from the Scani-valve and the PPB systems were reduced and tabulated on-line 
utilizing the AEDC-PWT digital computer facilities. Also, the static and total pres- 
sure parameters were stored together with identification and test condition parameters, 
on 800 bpi, IBM compatible digital magnetic tape for post-test analyses. 

The microphones (or pressure transducers) used for this test were Kistler 601 L quartz 
pressure transducers. These microphones have a resonant frequency of 130 kHz and 
previous calibrations conducted at NASA-MSFC indicate that frequency and phase 
were linear up to 20 kHz. The full scale pressure range was 300 psi and the active 
area of the diaphragm was 5/32-inch diameter. Kistler 553A charge amplifiers were 
used to convert the quartz charge to a voltage signal . The charge amplifiers limited 
the low frequency response to approximately 10 Hz. The amplifiers were also housed 
inside the model. A typical microphone installation is shown in Figure lib. 

The Kistler pressure transducer systems were calibrated in place by the use of a 
Photocon model PC-125 acoustic calibrator which was modified so as to be hand 
held over each pressure transducer. The calibration level was 150 dB (referenced 
to 0.0002 dynes per square centimeter) at 1 kHz. 

The data acquisition system for the fluctuating pressure measurements consisted of the 
NASA-MSFC, 20 kHz constant bandwidth multiplexer and 14 track, 1 .5 MHz 
magnetic tape recorder. A block diagram of the system is shown in Figure 12. The 
multiplexer provided close phase matching for data groups recorded on a single 
track and also enabled up to 122 data channels to be stacked on fourteen tracks. 

The data channels were stacked 7 and 9 to a tape track with each channel having 
a 10 Hz to 20 kHz frequency response. To provide for cross-spectrum analyses between 
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data recorded on different tape tracks, it was necessary to record a crystal oscillator 
signal and a wide-band white noise signal on all channels at the beginning of each 
reel of tape since phase matching across tape tracks was limited by the alignment of 
the tape recorder heads and type skew. Also, the output from the crystal oscillator 
could be used in conjunction with a reference modulator to correct for tape speed 
errors. These signals were recorded continuously with the microphone data during 
the test. A complete system calibration was conducted immediately prior to the test 
and spot checked periodically during the test. 

The vertical position of the various protuberances was measured by a linear poten- 
tiometer which was housed inside the model . The output from the potentiometer was 
conditioned through the AEDC-PWT force and moment readout system (FAMROS) 
and subsequently reduced and tabulated on-line. Calibration of this instrument was 
performed prior to the start of the test and checked periodically during the test. 
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3.0 


TEST DESCRIPTION 


3.1 Test Conditions 


Basic tests of all protuberance configurations were conducted in the Mach number 
range from 0.60 to 1 .60 and at a constant Reynolds number of 3 x 10 6 per foot. 

This unit Reynolds number was near the maximum that could be held constant in 
Tunnel 16T for the specified Mach number range. The 4-inch diameter protuberance 
was tested also at a Reynolds number of 1 .5 x 10 6 per foot and the 8-inch diameter 
protuberance at Reynolds number of 1 .5 x 10 6 and 4.5 x 10 6 per foot. For the 4.5 
million Reynolds number condition, the maximum Mach number attainable was 1 *2; 
however, data were obtained at Mach number 1 .4 at a slightly reduced Reynolds 
number. Variations of unit Reynolds number with Mach number are presented in 
Figure 13a for the various configurations. Similar variations of dynamic pressure, 

^oo' W '^ ^ach number are presented in Figure 13b. 

Tests were conducted for the 4-inch diameter protuberance at a unit Reynolds 
number of 3 x 10 6 per foot with an artificially thickened boundary layer over the 
test panel . For this study, the boundary layer was tripped near the leading edge 
of the model so that the effects of boundary layer thickness on the protuberance flow 
field could be evaluated. 

During the latter part of the test, experiments were conducted to develop effective 
flow visualization techniques in the region affected by the protuberance flow field. 

These experiments were performed concurrent with other phases of the test by using 
a 2-inch diameter, fixed protuberance located near the edge of the model, off the 
test panel . The technique evolved in these experiments was then used to obtain 
photographs of the flow patterns around the 4-inch diameter protuberance installed 
on the test panel . For the test, the Mach number ranged from 0.60 to 1 .60 with 
unit Reynolds number held constant at 3 x 10 6 per foot. For all these studies, oil 
with paint pigment mixtures of various colors was released forward of, around, and 
behind the protuberance through small orifices so that the various regions of the 
perturbed flow field could be made visible. The individual colors of oil were con- 
trolled with solenoid valves which were regulated from the control room. The f 

reservoir containers for the oil were pressurized to facilitate the flow of oil over 
the test panel surface. The resulting flow patterns were recorded on color still and 
motion pictures using cameras positioned on the top wall of the tunnel test section 
at a location approximately over the protuberance. The primary oil flow studies, 
using the 4-inch diameter protuberance, were conducted during the last two days of 
testing. For these studies, the microphones were removed from the model, and the 
holes were sealed with epoxy. Further, the test panel was painted with white enamel 
to improve color identification of the oil flow patterns. This technique proved 
extremely efficient for flow visualization over the test panel since it eliminated the 
need to interrupt testing between test conditions as required for conventional flow 
visualization studies. That is, the oil for the present technique remains sufficiently 
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fluidic that it adjusts to a change in flow pattern rather easily; whereas, conventional 
oil flow techniques (titanium dioxide solutions for example) require that the oil 
evaporate leaving a residue of substance in order to identify the flow pattern. A 
titanium dioxide solution was applied to the test article for the last run of the test 
to identify the flow pattern over the wall of the 4-inch diameter protuberance at 
Mqq =1.0, h/D = 2.0. For this run, the surfaces of the test panel and protuberance 
were painted black. 

Test Procedure 

Prior to testing on each day, a complete system calibration was performed. This 
"pre-test" calibration consisted of zeroing the protuberance height potentiometer 
readout, leak checking of all the pressure orifices and probes, and both acoustically 
and electrically calibrating the microphone/ data acquisition system. The acoustic 
calibration consisted of exciting each microphone at 150 dB with the Photocon PT-125 
calibrator, and adjusting the input electronics to the tape recorder so that equal 
sensitivity was obtained for each microphone channel . The electrical calibration 
consisted of inserting a 458 mv (equivalent to 170 dB which was full scale) 1000 Hz 
sinusoidal oscillator signal into the microphone charge amplifier self-test inputs and 
recording the resulting signal on the data tape. This calibration was followed by a 
45.8 mv (equivalent to 150 dB) 1000 Hz self-test input which was recorded also on 
the data tape. Further, a broad band random noise signal was recorded on tape 
through the self-test input. The sinusoidal self-test calibration was repeated at the 
beginning of each new reel of data tape during each day of testing. Periodically, 
the entire microphone/data acquisition cdlibration sequence was repeated following 
a day of testing for hot "post-test" calibration records. 

For air-on testing of the cylindrical protuberances, the protuberance height was 
varied with Mach number held constant for a given run. In general, data were 
obtained for monotonic increases in protuberance height because of possible hysteresis 
effects associated with boundary layer separation. For several configurations and test 
conditions, data were obtained aiso for variations in protuberance height which 
monotonically decreased so that hysteresis effects could be evaluated. There was no 
noticeable effect on the test results. The procedure of varying protuberance height 
while holding Mach number constant was adopted because of the lesser time involved 
in protuberance height variations as compared to variations in Mach number. Howeve 
for the RCS and APS protuberances, variations in height were not required, and data 
were obtained for monotonic variations in Mach number - with both increasing and 
decreasing trends. 

The test procedure for each cylindrical protuberance was as follows: 

• Data acquisition with the protuberance and forward rake in 
the fully retracted, flush, position. 


• Data acquisition with the protuberance in the retracted 
position and the forward rake in the fully extended postion. 

• Data acquisition with the rake in the retracted position and 
the protuberance extended to heights ranging from zero to 
eight inches in fixed increments of 0.25 diameters. 

The data acquisition procedure consisted of, first, acquisition of the static pressure 
data followed by acquisition of the fluctuating pressure recordings. Previous experi- 
ence with the present instrumentation systems indicated that noticeable instrumentation 
noise was generated by the operation of the Scani-valve cyclic drive systems which 
affected the microphone signals. Thus, static pressure data and microphone data were 
obtained in sequential order rather than simultaneously. This procedure also allowed 
for examination of the on-line tabulation of mean flow data prior to proceeding to 
new test conditions. 

For several configurations and tunnel test conditions, data were also obtained for 
transient variations in protuberance height (with Mach number held constant) and 
transient variations in Mach number (with protuberance height held constant) . For 
these runs, the protuberance height and Mach number were monitored and static pres- 
sure data were obtained at discrete intervals without stopping the transient run for 
data acquisition. The microphone outputs were recorded continuously for the duration 
of each transient run. 
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4.0 

4.1 


4.2 


4.3 


DATA REDUCTION PROCEDURE 


Introduction 


It was necessary to perform a relatively extensive reduction of the measured data to 
define the tunnel test environment and the basic mean-flow and unsteady-flow 
environments induced by three-dimensional protuberances at transonic speeds. It is 
convenient to discuss the data reduction in terms of 1) tunnel flow parameters, 

2) mean-flow parameters for the protuberance flow field, and 3) unsteady-flow 
parameters for the protuberance flow field. It should be noted that the computation 
of mean flow parameters was performed on-line during the test; whereas, the compu- 
tation of the unsteady-flow parameters was performed off-line, subsequent to the test, 
utilizing the microphone recordings stored on magnetic tape. 

Tunnel Flow Parameters 

The computation of tunnel flow parameters was performed on-line during the test 
using standard data acquisition and computer program routines adopted for this 
purpose by the engineering personnel at AEDC. It suffices to state that the following 
flow parameters were computed and tabulated utilizing the facilities at AEDC. 


q a» 

P co 
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• Free-Stream Mach Number, 

• Free-Stream Dynamic Pressure, 

• Free-Stream Static Pressure, 

• Free-Stream Unit Reynolds Number, 

• Stilling Chamber Total Pressure, 

• Stilling Chamber Total Temperature, 

Mean- Flow Parameters 


The computation of mean-flow parameters over the test article for the various 
configurations consisted, primarily, of reducing the measured static and total 
pressures to descriptive engineering terms. The static pressure measurements 
were reduced to ratios of the local to free-stream pressure, P^/P^ and to 

standard pressure coefficients, C . The total pressures from the boundary layer 

P 

rakes were reduced to ratios of velocity, U n /U Q0 , and to profile Mach number, M n . 

The boundary layer displacement thickness, 5*, and momentum thickness, 0, were 
computed also from the total pressures at both the forward and aft rake positions. 
Further, the boundary layer displacement thickness and momentum thickness were 
computed at model station zero utilizing Bie's empirical relationships (Reference 22) . 
The reduction of all static and total pressure measurements to mean-flow parameters 
were performed utilizing the digital computer facilities in the AEDC Propulsion Wind 
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Tunnel Facility. Pertinent equations used in the computation of the mean-flow 
parameters are presented in Appendix A. A typical printout of the on-line data 
is presented also in Appendix A. 

4.4 Unsteady Flow Parameters 

The computation of unsteady flow parameters to define the salient fluctuating pres- 
sure characteristics of the perturbed flow field induced by three-dimensional protu- 
berances consisted of reducing the microphone recordings to descriptive statistical 
parameters. The procedure for this data reduction was very complex and certainly 
the most time consuming phase of the study. Both analog and digital techniques were 
employed in the data reduction process as noted in the following sections. 

4.4.1 Analog Data Reduction 

The analog data reduction instrumentation consisted of the following: 

• A 1 .5 MHz magnetic tape recorder with a nine channel, 
demodulator/tuner system for demultiplexing the FM signals 
recorded on each tape track. The demultiplexed signals 
were a reconstruction of the originally recorded microphone 
signals for the frequency range from 0 to 20 kHz. 

• A Bruel and Kjaer Type 2112 Audio Frequency Spectrometer. 

The B&K spectrometer provided for 1/3-octave spectra and 
overall sound pressure level analyses. 

• A Bruel and Kjaer Type 2305 Graphic Level Recorder. The 
B&K Level Recorder provided plots of the 1/3-octave spectra 
and overall sound pressure level. 

The analog instrumentation was used in the reduction of all fluctuating pressure 
measurements to overall sound pressure level (OASPL) which were plotted on a dB 
scale by the graphic level recorder. Approximately 15 seconds of the real-time 
microphone signals were used in the analysis of each OASPL. Further, the frequency 
range covered by the analog equipment for the determination of the OASPL was from 
2 Hz to 20 kHz. 

The procedure adopted for the OASPL data reduction was to, first, adjust the data 
reduction electronics to give a 170 dB reading on the graphic level recorder for a 
458 mv calibration signal recorded on each reel of tape for each data channel to be 
reduced. The calibration signal was recorded for each channel so that a permanent 
record could be maintained. Following the calibration set-up, the desired data 
channels were reduced to OASPL. This procedure was repeated for data on each new 
reel of tape. The OASPL data were reduced further to root-mean-square coefficients 
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of fluctuating pressure, AC (RMS), by hand computation, using the following 

r 

equation to relate the OASPL to the RMS fluctuating pressure: 

OASPL =20 log 10 AP(RMS) + 127.582 dB (4.1) 

where the OASPL is referenced to 2 x 10“ 5 Newtons per square meter and 
AP(RMS) is in units of pounds per square foot. 

The analog instrumentation facilitated a rapid reduction of the fluctuating pressure 
measurements to RMS coefficients, thus, enabling the important microphone channels 
to be selected for digital reduction to more complex statistical parameters. The analog 
data reduction instrumentation was employed also to obtain 1/3-octave spectra of 
representative microphone channels for comparison with digital computations. 

4.4.2 Digital Data Reduction 

Selected data were reduced to descriptive terms using a CDC 3300 digital computer 
facility. The output from the magnetic tape system was converted to digital form 
using the CDC 3300 computer equipped with Texas Instruments Model 846 analog-to- 
digital converters and two model 845E multiplexers. The computer system was capable 
of acquiring analog data and converting it to digital format over two simultaneous 
channels at a rate of 41 ,700, 12-bit samples per second per channel . The converted 
data was recorded on 7-track digital magnetic tape for subsequent analysis and tabu- 
lation. A comprehensive description of the computer programs used for digital analysis 
of the fluctuating pressure data is presented in References 23 and 24. The computations 
consisted of the following parameters. 

• Root-Mean-Square Fluctuating Pressure Level, AP(RMS), psf 

• Overall Sound Pressure Level, OASPL, dB 

• 1/3-Octave Center Frequency, f, Hz 

• 1/3-Octave Reduced Frequency Parameter, f£/U c 

• 1 /3-Octave Sound Pressure Level, S(f), dB 

• Power Spectral Density at 1/3-Octave Center Frequencies, <|>(f), (psf) /Hz 

• Normalized Co-Spectral Density, c(£,f) 

• Normalized Quad-Spectral Density, q(£,f) 

• Total Cross- Spectra I Density, G(£,f), (psf) 2 /Hz 

• Phase Angle, a!(£,f), deg 

• Coherence Function, y (£,f) 

• Convection Velocity, U c (f), ft/sec 


14 


The analog forms of pertinent data reduction equations for the foregoing unsteady 
flow parameters are presented in Appendix B. 


It was noted in the preceding section that the fluctuating pressures were reduced to 
RMS levels using analog instrumentation. However, data which were selected for 
more complex analysis on the computer were also reduced to RMS levels, thus, pro- 
viding a means of cross-checking the results. A discussion of the comparison between 
analog and digital results is presented in Section 5.3.3. The digital computation of 
RMS fluctuating pressure was performed using digital samples taken at N points in 
time as follows: 


AP(RMS) = 


N 

- E 
N & 


n 


£ 


(4.2) 


where P n is the digital sample of the pressure P(t) with P = 0. 


The method of digital computation of the spectral density parameters was based on a 
direct, finite, Fourier transform of the pressure time signals which were thence con- 
verted to power spectral density. The "Fast Fourier Transform" (FFT) algorithm was 
used to calculate the finite Fourier transform. The method is briefly summarized in 
the following paragraphs. 


Let X(j) and Y(j) be the microphone pressure signals recorded at times j At, that is, 
the time functions defined at N points in time. At apart. The finite Fourier transform 
of X and Y are: 


A(k) 


JL 

N 


N - 1 

2 X(j) exp y- 2ir i 


(4.3) 


j = 0 


B(k) = 


N 


N - I 

Y (j) ex P 

j~0 


(* 2,i TJ ) 


(4.4) 


where 


k = 0,1,2, ..., N-l 


i = 




A spectral window is required and may be applied by weighing the X and Y samples; 
or, for the window used in the present analysis (see Reference 24), may be applied by 
smoothing the real and imaginary parts of A and B with weights -1/4, 1/2, -1/4. Thus, 
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modified periodograms are used (a periodogram obtained from data to which a window 
other than the Dirac comb function has been applied is known as a modified periodo- 
gram) . The modified direct- and cross-peri odograms of the signal pairs are: 



I x (k) = | A’(k)| 2 

(4.5) 


Iy(k) = I B'(k) I 2 

(4.6) 


I xy (k> = A'(k) B' * (k) 

(4.7) 

where A'(k) and B'(k) are the smoothed forms of A and B and the asterisk means 
complex conjugate. These periodograms may be converted to power spectral density 
with each value having 2 degrees of freedom. To increase the statistical accuracy, 
a series of p periodograms from adjacent sets of data are averaged, thus giving power 
spectral density results with 2p degrees of freedom. The method is to divide the N 
values of X(j) and Y(j) into p sections each of length m, sum the I(k) for each k and 
divide by p, i .e . , 


P 

T x (k) =4 Z [i x (k)l 

(4.8) 


P 

yk)=^E n y (k)i r 

H r = 1 

(4.9) 


P 

Ixy(k) = ~ Zi f*xy(k)] r 
P r = l 

(4.10) 

where now k = 0,1,2, . . ., m - 1 and the total number of values of each signal which 
has been used is mp. Power spectral density results are defined for m/2 + 1 positive 
frequencies from 0 Hz to 1/(2 At) Hz, giving a frequency interval of l/(m At) Hz. 
The direct and cross-power spectral densities are given by 


<p (k) = 32 m At I x (k) 

(4.11) 


<|> y (k) = 32 m AtT y (k) 

(4.12) 


<P xy (k) = 32 m AtT xy (k) 

(4.13) 


for the spectral window defined for the modified periodogram as discussed in 
Reference 24. 
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The power spectral density expressions given above were obtained at the frequency 
interval of l/(m At) Hz and subsequently converted to values representative of 1/3- 
octave bands. USAS! preferred frequencies were used in the specification of the 1/3- 
octave frequency intervals and center frequencies. This averaging process was per- 
formed utilizing power spectral densities from Equations 4.11 to 4.13 contained within 
each 1 /3-octave band as follows: 


♦ x « - A? / E * x ( k > Af k ■ *( f > 

/3 f t 
, fu 

y f > - at^t E y k > Af k - yo 

/3 f t 

f u 

Af k = G(|,f) 

/ 3 H 


(4.14) 


(4.15) 


(4.16) 



where Afy represents the 1/3-octave band, Af|< is the interval over which the 

power spectral densities values tj>(k) apply, and f^ and f are the lower and upper 
frequencies, respectively, of the 1/3-octave band, Afy . At the edges of each 

1/3-octave band, the value of $(k) was obtained by linear interpolation between the 
appropriate neighboring values of $(f). It should be noted that this averaging process 
results in average energy per Hz within each 1/3-octave band rather than the total 
energy in the 1/3-octave band. Further, this averaging process increases the number 
of degrees of freedom of the direct and cross-power spectral density estimates in the 
same way as averaging the periodograms as indicated in Equations 4.8 to 4.10. See 
Section 5.3.3 for a discussion of the statistical accuracy of the computation. It will 
be noted in Section 5.3.3 that, to further increase the statistical accuracy of the 
computations at the low frequencies (below 200 Hz) and concurrently maintain 
acceptable computational times, it was necessary to average power spectral density 
values over full octave bands. Thus, below 200 Hz, data are presented only at the 
full-octave center frequencies. 

The method of calculating direct and cross-power spectral density results in which 
the frequency range is divided into a series of bands, as described in Reference 24, 
was used in the data reduction computer programs. This method requires the use of a 
digital filter to successfully halve the frequency content of the time signal . The 
iurgest frequency interval is obtained in the upper frequency band and the interval is 


I 
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successively halved in each lower band. This process accounts for the nonlinear 
variations of the number of statistical degrees of freedom and normalized standard 
error with frequency as discussed in Section 5.3. A considerable amount of computer 
time was saved by the digital filtering process compared with the method of calculating 
the power spectral density with a constant frequency interval low enough to obtain the 
one-third octave band centered at 25.1 Hz. 

The direct and cross-power spectral densities were used in the calculation of several 
useful functions which express the properties of the fluctuating pressure. Noting 
that the cross-power spectral density consists of both real and imaginary components, 
such that 

G(£,f) = C(f ,f) - i Q(£/f) (4.17) 

the following functions can be defined: 

Normalized Co-Spectral Density 

C(f ,0 

c(|,f) = — r (4.18) 

rtff)^)]* 

where C(£ ,f) is the real component of the cross-power spectral density, G(£ ,f). 


Normalized Quad-Spectral Density 


q(£,f) 


Q(M) 


(4.19) 


where Q(£,f) is the imaginary component of the cross-power spectral density, G(£,f). 
Coherence Function 


Phase Angle 


r 2 («/f) 


|G(l,f)| 2 



(4.20) 


(4.21) 
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Narrow-Band Convection Velocity 


u c (f) 


2tt f £ 

a(£,f) 


(4.22) 


Broad-Band Convection Velocity 


k 

£ U c(f) n Af„ 

U c = (4.23) 

£ Af " 

n = 1 

A typical printout of the digital computations is presented in Appendix B. Also, the 
analog form of pertinent equations used in the computations are discussed in 
Appendix B. 
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5.0 


PRECISION OF. MEASUREMENTS 


5.1 Tunnel Flow Parameters 

The uncertainties associated with maintaining tunnel conditions were computed by 
AEDC personnel (Reference 25). These uncertainties are estimated as follows: 


Mach Number, Subsonic 

±0.005 

Mach Number, Supersonic 

±0.010 

Total Temperature 

±5°F 

Total Pressure 

± 5 psf 


The uncertainties in tunnel flow parameters result from variations in setting and 
maintaining tunnel conditions. The longitudinal variation of Mach number along the 
centerline of the test section in the vicinity of the model is not included in the above 
values and reaches a maximum of ± 0.007 at supersonic Mach numbers. 

5.2 Mean-Flow Parameters 


The pressure measurements of the model pneumatic system have been evaluated by 
measuring a common pressure on all pressure scanning switches (Scani-valves) and 
on the tunnel Precision Pressure Balance (PPB) system. The variations were a function 
of the test condition total pressure (or Unit Reynolds number) and Mach number. The 
variations were linear with tunnel total pressure such that, when expressed in terms of 
pressure coefficient, the variations were constant at a given Mach number. The maxi- 
mum variations in Cp, based on the worst cases encountered, ranged from ± 0.009 at 
M^ = 0.60 to ± 0.026 at M^ = 1 .60. These variations were computed directly from 

the variations in measured static pressures using the PPB system as a standard and do not 
take into account the uncertainties associated with either the PPB system or the 
uncertainties in dynamic pressure. 


In addition to the uncertainties presented above, variations in C were measured over 

P 

the test panel for protuberance-free flow conditions (protuberance in fully retracted 
position). These spatial variations resulted from nonuniform flow over the test panel 
and are attributed to the aerodynamic characteristics of the model rather than varia- 
tions in tunnel flow. A discussion of the spatial variations in Cp due to nonuniform 
aerodynamic effects is presented in Section 6.2.1 . 


5.3 Unsteady Flow Parameters 


Precision in the measurement and reduction of fluctuating pressure data is an 
important, but often neglected, characteristic of the results. Almost any review of 
published test results for fluctuating pressure environments will reveal considerable 
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scatter (and often a lack of agreement) between the results from the various experi- 
mental studies. One good example is the case of attached turbulent boundary layers. 

Here, numerous studies have been conducted in both wind tunnel facilities and in 
free flight; however, considerable scatter has been shown in the results from the 
various sources. Further, only in rare cases has the precision of measurements been 
presented with the data. Thus, it must be concluded that the errors and uncertainties 
contribute significantly to the lack of agreement between experiments. This is 
particularly true for fluctuating pressure environments of relative low overall level 
since numerous sources of error can be present in the measured data which are of the 
same order of magnitude as the fluctuating pressure of the phenomena under study. 

For the present test results, a comprehensive analysis of the precision of measurements 
has been made. For the purposes cf this analysis, the precision of measurements will 
be distinguished by 1) errors, and 2 ) uncertainties. Errors are defined as deviations 
from the true values which are deterministic in the sense that they can be assessed 
through proper calibration and analysis. Further, errors may be applied directly to the 
measured data to minimize the inaccuracies in the final results. Uncertainties, on the 
other hand, are defined as deviations in the measured data about the true value and 
are nondeterministic in the sense that they cannot be applied directly to the data to 
determine the true value. Uncertainties represent the confidence band or confidence 
limits of the test results. At this point it should be noted that the data presented 
herein have not been corrected for the estimated errors . However, the errors are 
presented in a form which is consistent with that of the results, and, thus, may be 
applied with relative ease to selected results. The extent of the present data reduc- 
tion and analysis (with the primary aim of determining the salient features of the pro- 
tuberance flow fields) precluded the correction of all the data for known errors, and 
thus the data is presented in uncorrected form with the errors specified in this section. 

Inaccuracies in the measurement and reduction of fluctuating pressure data result 
from a number of sources ranging from the inaccuracies associated with the data 
acquisition electronics to the statistical inaccuracies of the final computations. The 
purpose of this section is to note the predominant sources of inaccuracies and to assess their 
order of magnitude. It will be noted that errors presented herein could be interpreted 
as uncertainties in the data, since no effort has been made to adjust the test results for 
known error contributions. For example, reasonable estimates of instrumentation 
noise, tunnel background noise and system response are determinable; however, these 
contributions are relatively small compared to the level of the fluctuating pressure 
associated with the protuberance flow field. Thus, the contributions from these 
extraneous sources could be regarded as uncertainty errors in the precision of the 
measurements if so desired. On the other hand, the data may be refined in accuracy 
by applying the appropriate corrections presented herein. 
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The following tabulation consists of errors and uncertainties which are known to be 
prevalent in the data presented herein: 

• Extraneous Noise Errors 

— Instrumentation System Noise 
— Tunnel Background Noise 

• Frequency Response Errors 

— System Frequency Response 
— Finite Microphone Size 

• Data Reduction Errors and Uncertainties 

— Statistical 
— Technique 

In the following sections, the errors and uncertainties associated with each of the 
foregoing sources are discussed and assessed. For the purposes of this analysis, it is 
assumed that the various extraneous noise sources are uncorrelated between microphone 
channels such that noise correlation errors are not present in the cross- spectra I density 
computations. A thorough examination of the noise correlation problem is in progress; 
however, these results are not available at this time for inclusion in the present report. 
It is anticipated that a separate report will be prepared in the near future defining in 
more detail the uncertainties in the cross-spectra due to the possible correlation of 
extraneous noise. The present precision of measurement analysis has considered the 
overall fluctuating pressure level and the power spectra. This analysis is summarized 
in Section 5.3.4. 

5.3.1 Extraneous Noise Errors 

The presence of extraneous noise in recorded data may have a significant effect on 
the accuracy of the various statistical parameters used to describe the characteristics 
of random phenomena under study. For the present wind tunnel investigation, the 
recorded fluctuating pressure data contained extraneous signals superimposed on the 
true pressure-time histories. The notable extraneous noise sources are the instrumen- 
tation system noise associated with the data acquisition equipment, and aerodynamic 
noise associated with the test facility air stream. Since the frequency domain of 
interest for the true fluctuating pressure data may coincide with that of the extraneous 
noise, removal of the noise cannot be accomplished by a simple filtering of the 
measured signals. Thus, one must resort to statistical means of extracting the true 
fluctuating pressure data from the measured data* One other alternative is to evaluate 
the statistical characteristics of the extraneous noise and include these data as errors 
in the precision of measurements. This latter approach is certainly feasible when the 
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fluctuating pressures associated with the unsteady flow phenomena under study are 
much larger than the effective level of the extraneous noise. For the present test/ 
the levels of extraneous noise are relatively small compared to the fluctuating pres- 
sure levels associated with the protuberance flow field and, consequently, the 
extraneous noise may be considered as uncertainties in the precision of measurements. 

A brief assessment of the two predominant extraneous noise sources - instrumentation 
system noise and tunnel background noise - are presented herein. 

Instru me ntation Noise 

The instrumentation system noise was analyzed from the wind-off data records. The 
overall system noise levels expressed as equivalent RMS fluctuating pressure and 
representing the mean of the track-to-track deviations are presented as variations 
with VCO carrier frequency in Figure 14a. The presentation of results representing 
the mean of the 14 tc r 'C tracks is felt to be appropriate here, since track-to-track 
deviations about the mean are effectively errors in the error and, thus, have a second 
order effect on the precision of the measurements. The system noise increased with 
increasing carrier frequency. This increase is explained by: 1) the noise of the tape 
recorder increases at higher frequencies corresponding to the higher carrier frequencies, 
2) the high frequency channels are more susceptible to tape flutter due to their small 
percentage deviation, and 3) the harmonics of the lower frequency channels are 
present in the higher channel pass bands. Previous calibrations conducted by NASA- 
MSFC indicated that worse case signal-to-noise ratios run from 26 dB to 38 dB relative 
to the full scale deviation of all channels. These data are consistent with the present 
test results which show mean signal-to-noise ratios ranging from 27.7 dB for VCO 
channel No. 9 to 37.3 dB for VCO channel No. 1 . 

The power spectra of the instrumentation noise for the various microphone channels 
examined were all similar in shape. Thus, it suffices to examine only the mean power 
spectrum obtained by averaging the spectra for the nine VCO channels. The tape 
tracks from which the nine VCO channels were selected were typical of all 14 tracks. 

It should be noted that variations in instrumentation noise from channel to channel 
results in a corresponding variation in measured data. Thus, c comparison of averaged 
instrumentation noise with averaged measured data introduces little error into the final 
analyses. Further, deviations about the average due to channel-to-channel and track- 
to-track variations are errors in rhe error and thus mcy be considered as having a 
second order effect as previously mentioned. The mean spectrum for the instrumentation 
noise is presented in Figure 14b as the distribution of equivalent mean-square fluctuating 
pressure with frequency. 

Tunnel Background Noise 

The most difficult extraneous noise source to evaluate is the wind tunnel background 
noise. Accurate measurement of the wind tunnel background noise is complicated by 
the fact that any calibration probe will generate its own aerodynamic flow field which 
may result in the measurement of erroneous tunnel noise. In the past, little attention 
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has been given to the problem of wind tunnel noise, and the general trend has been 
to ignore the problem. Fortunately, measurements have been made for the test 
facility used in the present study; although, these data do not completely define the 
noise environment. The data, presented in Reference 21, consist of overall RMS 
fluctuating pressure and power spectra,. These data were obtained using a slender, 
modified cone having a 5-degree half angle. Microphones were flush mounted along 
one meridian of the cone and spaced so that measurements could be obtained in both 
the laminar and turcnlent regions of the cone boundary layer. Typical overall 
fluctuating pressure levels, in terms of the ACp(RMS), are presented in Figure 15a 
for the Mach number range from 0.60 to 1 .40. Maximum free-stream fluctuating 
pressures were present at Mach numbers near 0.70. It will be noted that data for the 
present study as obtained from the test panel with the protuberance in the fully 
retracted position are presented also in Figure 15a, and these data compare reasonably 
well with the calibration results. The data for the test panel at each Mach number 
represent the mean of approximately 88 microphone readings taken along the center- 
line of the test panel . As for the instrument noise case, the deviations about the 
mean have a second order effect and, thus, are neglected. The agreement between 
the calibration probe results and the test panel results is not too surprising even though 
a relatively thick turbulent boundary layer was present during the test panel measure- 
ments. It will be shown that the fluctuating pressures contributed by the attached 
turbulent boundary layer range from 10 to 20 percent of the total fluctuating pressure 
level based on a mean square comparison. Based on this finding, it is concluded that 
the protuberance-free, test panel results can be used to estimate the extraneous noise 
by applying an empirical correction to these data for the turbulent boundary layer 
noise. That is, with the protuberance in the retracted position (denoted by the sub- 
script, 0) consider the measured pressure, P(t) 0 , as the sum of three stationary random 

processes - tunnel noise, P(t)y, instrumentation noise, P(t)j, and attached boundary 
layer noise, P(t)g, where the total extraneous noise, P(t)^, is the sum of P(t)y and 
P(t)j, such that 

P(t)o = P(t>T + P «I + P«B = P «E + P «B (5-1) 

It follows that, if the various components of the measured signal, P(t) 0 , are 
statistically independent, then the mean square fluctuating pressure is given by 

[P(RMS)] 2 = e f AP(RMS)1 ^ + e[AP(RMS)] 2 + [AP(RMS)}| 

= e [AP(RMS)I 2 + [AP(RMS)] 2 (5.2) 

where the prefix e denotes an error term. 
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In RMS coefficient form, the mean square fluctuating pressure as given in Equation 
5 .2 becomes 

1 

[ ACp(RMS)] 0 = |e [AC p (RMS)] 2 T + e [AC p (RMS)], + [AC p (RMS)l 2 J 2 

1 

= {e [AC p (RMS)] 2 E + [AC p (RM5)] 2 B j- 2 (5.3) 

Thus, the test panel results may be used to obtain an estimate of the total extraneous 
noise level by direct subtraction of the mean square boundary layer noise from the 
mean square test panel data. A comprehensive analysis of turbulent boundary layer 
fluctuating pressures has been made by Lowson (Reference 26). This analysis was 
based on a broad selection of experimental results from a number of different sources. 
An empirical formula for the boundary layer fluctuating pressure level which appears 
to agree with the general trend in the data is given by 

[ ACp(RMS)] B = 0.006/(1 +0.14 M 2 ) (5.4) 

The variation of the estimated boundary layer noise from Equation 5.4 with Mach 
number is presented also in Figure 15a. 

Utilizing Equations 5.3 and 5.4, the following expressions are derived for the total 
extraneous noise and the tunnel noise. 

r ( 2 3.6 x 10 -5 ) 5 

e [ AC (RMS) ]_ = | f AC (RMS) ] 7 > (5.5) 

P E ( P (1+0.14M 2 ) 2 ) 


( 2 2 3.6 x 10“ 5 ) 2 

e ( AC (RMS) ] T = { [AC (RMS) ] - e [AC (RMS)]' r > (5.6) 

H K Pi (1 +0.14 M 2 ) ) 

where the instrument noise, e [AC (RMS)]t, may be obtained from Figure 15b. 

pi 

Variations of the total extraneous noise, instrumentation noise, and tunnel noise with 
Mach number are presented in Figure 15b showing a comparison with the total measured 
results. It is evident that, for the protuberance-free case, the tunnel noise dominates 
the measured signal from the relative magnitude of the contributing processes. It is 
also evident that little insight into the attached turbulent boundary layer case could 
be gained from the test panel measurements with the protuberance retracted. Again, 
it should be emphasized that the primary objective of this study was the examination 
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of the protuberance flow fields and not the undisturbed boundary layer noise. For 
this reason, the protuberance-free data are presented only as an indication of the 
extraneous noise levels in the protuberance flow fields in subsequent figures. 


The errors in the power spectra resulting from extraneous noise may be analyzed in 
a manner similar to that used in the overall fluctuating pressure level analysis. 
Again, it is necessary to rely on an empirical formula for adjusting the measured 
test panel results for the contribution of attached turbulent boundary layer noise. 
For the boundary layer, Lowson (Reference 26) derived the following empirical 
formula for the power spectral density. 


[q>(f)] B = [AC p (RMS)] 2 B 



2it f 




(5.7) 


where 



It will be noted that the power spectral densities for turbulent boundary layers are 
a function of free-stream dynamic pressure (q^), free-stream velocity (U^), and 

boundary layer thickness (5). There is much uncertainty regarding the appropriate 
boundary layer parameter to be used in the normalization of the power spectral 
density. Lowson used the thickness (8); whereas, others have used displacement 
thickness (8*), momentum thickness (9), and wall shear stress (t) . For the present 
report, 8* has been selected so that the results would be consistent with the two- 
dimensional protuberance results of Chyu and Hanly (Reference 27). However, it 
should be noted that protuberance diameter definitely controls the scale of the 
perturbed How field in the three-dimensional case, and any generalization of the 
data should include this parameter. 


For the analysis presented herein, comparison between two- and three-dimensional 
protuberance flows are discussed, and 8* appeared to be an adequate normalizing 
parameter. For the undisturbed flow case, it is felt that Equation 5.7 is appropriate 
for estimation of the boundary layer noise spectra even though 8* is used as the 
normalizing parameter. 


Using arguments similar to that for the overall fluctuating pressure analysis, the 
measured power spectral density for the protuberance-free test pane! may be expressed 
as: 
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r<Xf)J 0 = e [$(f)] T + e [t(f)]j + tt(f)] B 


= e[$(f)] E + [<|>(f)] B 


(5.8) 


where, as before, the subscripts T, I, B and E represent the tunnel, instrumentation, 
boundary layer and total extraneous noise contributions to the measured signal (sub- 
script 0). The errors in the measured power spectra are given, for the total extraneous 
noise and the portion due to tunnel background noise, by: 

e [<|>(f)] E = [*(f)] 0 “ W)] B 

e [<f»(f)l T = fcKf)] 0 - [<K01 b " e ft(f)3j (5.9) 


where [ <p(f ) ] D is defined by Equation 5.7 and e [<{>(f)] is given in Figure 14b. Power 
B 1 


spectra, representing the protuberance-free measured data, total extraneous noise, 
instrumentation noise, tunnel noise and boundary layer noise, are presented in 
Figures 16a through 16f for free-stream Mach numbers of 0.60 through 1 .60, respec- 
tively. It is noted in the spectra that the difference between the extraneous noise 
and the protuberance-free measurements cannot be distinguished since the contribution 
of the boundary layer noise is negligibly small. Further, the extraneous noise is 
attributed primarily to the tunnel background noise as was noted in the overall level 
results. 


At this point it should be noted that the tunnel background noise was affected also by 
the frequency response characteristics of microphone system (see next section). How- 
ever, it is convenient to leave the tunnel noise error (and total extraneous noise error) 
in a form which is uncorrected for frequency response errors as discussed in the next 
section . 

5.3.2 Frequency Response Errors 

In addition to the extraneous noise errors discussed in the preceding section, another 
important uurce of error should be evaluated - the frequency response. The frequency 
response of the data acquisition system is a function of 1) the electronics of the system, 
and 2) the attenuation characteristics of the microphone due to the finite size of the 
microphone sensitive area. In the evaluation of the frequency response errors, it is 
convenient to discuss, first, the effect of frequency response on the power spectra, 
and from this evaluation, assess the effect of the response error on the overall 
fluctuating pressure level . It will be noted that the frequency response errors cannot 
be expressed as an absolute value as was the case for the extraneous noise, but rather, 
must be expressed as a function of the measured power spectral density. Thus, each 
power spectrum presented herein will have an associate frequency response correction 
spectrum . 
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System Electronics 


The frequency response characteristics of the electronics of the microphone/data 
acquisition system are illustrated by results obtained from NASA-MSFC in Figure 17. 
The recording system frequency response was fiat for the frequency range from 20 Hz 
to 20 kHz. Thus, no errors are attributed to the recording system. The frequency 
response of the microphone electronics was obtained from discrete frequency calibra- 
tion signals inserted into the charge amplifier self-test inputs, recorded on tape, and 
subsequently played back for analysis. Typical calibration results are expressed as 
the ratio of the output-to-input power. The frequency response is shown to be flat 
within ± 1 dB for the frequency range from 20 Hz to 10 kHz. The frequency range 
of the present power spectra analyses corresponded to the 1/3 octave bands with 
center frequencies ranging from 25 Hz to 16,000 Hz, and the response attenuation 
was less than 3 dB for this frequency range. Before proceeding to the development of 
the error due to frequency response characteristics of the microphone electronics, it is 
convenient to assess the frequency response error attributed to the finite size of the 
microphone . 

Microphone Size Effect 


The measurement of a random pressure field by a microphone having a finite size 
sensitive area is subject to error because the microphone senses the spatial average 
over its face rather than a point value. To assess this error, it is necessary to have 
a priori knowledge of the true pressure field. Since the true pressure field is unknown, 
and in fact is precisely the parameter under study, determination of the microphone 
size effect becomes very difficult to assess. Several investigators have evaluated this 
problem through various mathematical models of the pressure field for attached boundary 
layers (References 28 - 33). However, their results are not directly applicable to the 
separated and wake flow fields as encountered in protuberance flow fields. In par- 
ticular, the correction factor for the finite microphone size is a function of the cross- 
spectrum characteristics of the pressure field. For the boundary layer case, sufficient 
data is available to provide an empirical formula for the cross-spectrum, and thus a 
priori knowledge sufficient to enable empirical size correction factors to be derived. 

The problem is considerably more complex for the three-dimensional protuberance 
flow field. First, previous data is not available with which to compare the present 
results. Second, if data were available, the nonhomogeneous character of the flow 
field would preclude the derivation of a general correction factor. Thus, only a 
crude approximation of the correction factor appears feasible without an extensive, 
and time consuming study of the problem. 

For the present analysis, it will be assumed that the turbulence structure for the pro- 
tuberance flow field bears some resemblance to that of the boundary layer. In par- 
ticular, it is assumed that the cross-spectra for the protuberance flow field has the 
same form as that for the boundary layer, and thus may be defined by the following 
empirical formula: 
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G(£,u) = 0(u) exp 


- C 


1 


0 ^1 


u 

' ' T 

c 

c 

- 


X exp 




(5.10) 


where ^ and £ 2 are the separation between points in the longitudinal and lateral 

directions, respectively, and C] and C 2 are constants which represent how fast the 
correlation decays in the longitudinal and lateral directions. For the boundary layer, 
Ci « 0.1 and C 2 ^ 0.7. For the protuberance flow field the decay appears to be 

more rapid in the longitudinal direction due to the dispersive nature of the flow field 
both upstream and downstream of the protuberance. In the neck region of the wake, 
the cross-spectra correspond closely to the exponentially decaying cosine with 
Ci > 0.2. In the separated region upstream of the protuberance, the decay is 

noticeably more rapid with the constant Ci ranging from approximately 0.4 to 0.8. 

It was noted in Reference 28 that an increase in the longitudinal decay rate would 
increase the sensitivity of the microphone, and consequently decrease the response 
error. Thus, a conservative estimate of the microphone size response error would 
correspond to the lower Cj value of 0.2 associated with the wake data. Lateral 
cross-spectra were not computed for the present data, however, it is reasonable to 
expect the lateral decay rate to be equivalent to the boundary layer (0.7) or less 
because of the three-dimensional character of the flow field. Transverse velocity 
components will have a tendency to decrease the lateral decay rate. Thus, a value 
of C 2 = 0.4 is arbitrarily chosen. For the range of Ci and C 2 which appears practical 
for the protuberance flow field, there does not appear to be a very significant change 
in the microphone correction factor (see Reference 28) . Thus, the selection of 
C] =0.2 and C 2 = 0.4 appears to be adequate for a prediction of the microphone cor- 
rection factor. The response function, M(f), which represents the ratio of the measured 
power spectral density to the true power spectral density in regard to microphone size 
effects, is presented versus normalized frequency in Figure 18. These data were 
obtained from Reference 28. Since the convection velocity is a function of frequency, 
the proper numerical values of U c may be used in estimating M(f) for a particular fre- 
quency. The radius of the microphone effective diameter, r, is 0.078 in. 

For a further approximation of the microphone size error, consider the neck region of 
the protuberance wake as in the preceding discussion. Here, the convection velocity 
(at high frequencies where microphone size becomes important) is approximately equal 
to the free-stream velocity. Thus, making this substitution for U c , the variation of 
M(f) with frequency may be approximated for each Mach number. These results are 
discussed in combination with the frequency response characteristics of the system 
electronics in the next section. 


Combined Response Error 


At this point, the frequency response functions, E(f) for the system electronics and 
M(f) for the finite microphone size effect, may be combined to give a total frequency 
response function, R(f) . Noting that: 

[Output Power] 

E(f) (5.11) 

[ Input Power] 

[Output Power] 

M(f) = — (5.12) 

[Input Power ] m 

and 

[Input Power ] g = [Output Power] (5.13) 

it foil ows that: 

[Output Power] 

R(f) = E(f) M(f) - (s.h) 

where [Output Power] is the power associated with the measured signal and 

[Input Power ] is the power associated with the fluctuating pressure field under 

examination. The total frequency response function, R(f), is presented in Figure 19 
and it is these data that should be used in the evaluation of the total frequency 
response errors. 

To utilize the frequency response function to evaluate the error which it contributes 
to the measured signal, it is necessary to re-examine in some detail the characteristics 
of the extraneous noise presented in the preceding section. The instrumentation system 
noise resulted primarily from the data recording electronics. Thus, the instrumentation 
noise was not affected by the frequency response characteri stics of the microphone 
system. It follows that the instrumentation noise error must be removed from the 
measured data prior to applying the frequency response correction. On the other hand, 
the tunnel noise error was affected by the frequency response of the microphone system, 
and, thus, it may be removed either before or after the frequency response correction. 
If the tunnel noise error is removed after the frequency response correction, then the 
tunnel noise correction must also be adjusted for the frequency response error. It is 
obviously more convenient to adjust the measured data for the tunnel noise error prior 
to applying the frequency response correction, and for this reason, the tunnel noise 
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(end total extraneous noise) were presented in the preceding section without the 
frequency response correction. 


Since the frequency response function is presented as the ratio of output-to-input 
power, the results may be applied directly to the power spectra adjusted for extra- 
neous noise errors. That is, the power spectral density, corrected for both frequency 
response and extraneous noise errors, is given by: 


0(f) = 


4tf)-e f0(f)! E 

R(f) 


(5.15) 


where e [ q>(t) ] ^ is given in Figure 16 and R(f) is given in Figure 19. Again it should 

be noted that R(f) as presented in Figure 19 is only an approximation due to the 
assumptions made regarding the cross-spectra of the fluctuating pressure field in the 
microphone size analysis. 


Equation 5.15 may be used to derive an expression for the total error in the measured 
power spectra due to the combined effect of the extraneous noise and the frequency 
response. Noting that 


e ft>(t)l E+R - <(>(0 ' 0(f) 


it follows from Equation 5.15 that 


e f 0(0 ] E _j_ R - <f>( f ) - 


{<Kf) - e [<Kf)l E } 


R(f) 


The error attributed only to the frequency '•esponse is given by 

ski 


1 1>(f) ] R = {l>(f)-e(t(f)) E } 


1 - 


(5.16) 


(5.17) 


(5.18) 


To assess the effect of the response errors on the overall sound pressure level, it is 
necessary to know the power spectra of the measured fluctuating pressure environment. 
Generally, the effect of the present response errors on the overall level will be 
negligibly small since the errors occur primarily at high frequencies where the power 
spectral density is small. However, noting that the overall fluctuating pressure level 
is related to the power spectral density by the following equation, the error may be 
computed for any given case. 
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AP(RMS) = 


1 

2 



#) df 


(5.19) 


For the present analysis, a broad band filter with a band pass from 2 Hz to 20 kHz was 
used in the reduction of the microphone measurements to overall RMS level . 


Thus, the error in the measured RMS fluctuating pressure, expressed in coefficient 
form, for the frequency range from 2 Hz to 20 kHz is given by 


e [AC p (RMS)] R 


f ACp(RMS)] 


20 kHz 



<K0-eItff)] E 

W) 



(5.20) 


where the subscript, R, denotes that the error is due to the system frequency response. 


It will be noted that the error resulting from the system response characteristics will 
always be negative, such that the measured AC p (RMS) is less than the true value. 

5.3.3 Data Reduction Errors and Uncertainties 

The primary errors and uncertainties in the final data result from inaccuracies 
associated with the data reduction process. Errors and uncertainties associated 
with the initial setup prior to acquiring the raw data are virtually eliminated since 
the calibration, as recorded on tape, are used in the data reduction process. Errors 
and uncertainties associated with the calibrations themselves are considered negligibly 
small as compared to those of the data reduction process. Thus, only the inaccuracies 
associated with the data reduction process are considered in the estimation of errors 
and uncertainties in the final data as presented herein. 

Data reduction errors and uncertainties may be classified as 1) statistical, and 
2) technique. The statistical inaccuracies result, primarily, from uncertainties 
associated with the computation of estimated values of the various statistical parame- 
ters. In random data analysis, use is made of finite bandwidths and finite sample 
lengths. Thus, the final computations are not true values, but estimations of the true 
values (see Reference 34). As bandwidth and/or sample length is increased, the 
uncertainties associated with the data reduction process are decreased. However, to 
obtain a certain degree of frequency resolution, finite bandwidths are employed. 

This, together with the practical necessity to limit the time of the sample, results in 
a certain degree of statistical uncertainty in the final statistical parameters. Tech- 
nique errors and uncertainties result from the inaccuracies associated with the equip- 
ment employed in the data reduction process and the ability of the person performing 
the data reduction. For example, the accuracy of setting up the data reduction 
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instrumentation, the accuracy of the various instruments employed, and the ability of 
the individual in reading the reduced data, say, from a graphic level recorder, will 
have a direct effect on the error and uncertainty of the final data. Technique 
inaccuracies are very difficult to assess, since several individuals, instruments, and 
sequences may be involved in the interim between the raw data as recorded on tape 
and the final data parameter as presented in the report. It should be noted that 
technique errors and uncertainties are not limited to microphone measurements; but, 
rather, are present in all experimental data. However, for the present precision of 
measurement analysis, only the microphone data will be considered in the assessment 
of technique errors and uncertainties. 

The purpose of this section is to define the inaccuracies associated with the computa- 
tion of each statistical parameter. In particular, the inaccuracies will be estimated 
for the RMS fluctuating pressure levels, the power spectral densities, and the cross 
power spectral densities. 

The analysis of the inaccuracies for these parameters are discussed separately in the 
following paragraphs. 

Overall Fluctuating Pressure Level 

The most convenient term to use in the specification of the statistical uncertainties in 
the overall fluctuating pressure level is the normalized standard error (e). The nor- 
malized standard error is the ratio of the standard deviation of the estimate to the true 
value. Physically, the normalized standard error may be interpreted as a fluctuation 
or scatter superimposed on the true estimate. In terms of confidence level, the 
probability that the computed (or estimated) value of AP(RMS) will be within ±e of 
the true value is 68 percent, or one has 68 percent confidence that the uncertainty 
will be within the limits imposed by ± e . The normalized standard error for the 
overall RMS fluctuating pressure level is given by (see Reference 34). 


e [ AP(RMS)] 

d 


s.d. f AP(RMS) I 
[ AP(RMS) ] true 



(5.21) 


where the prefix e implies the uncertainty normalized standard error of the term 
with the brackets, s.d. implies standard deviation; B is the bandwidth of the analysis, 
and T is the sample length of data used. For the present analysis of AP(RMS), using 
both analog and digital techniques, the BT product was sufficiently large that e was 
negligibly small . For the analog analysis: 


B « 20,000 Hz 
T = 1 5 sec 
BT = 300,000 

e [ AP(RMS)] $ * ! = 0.0009 (5.22) 

2 v'bT 

For the digital computation, the normalized standard error may he expressed directly 
as a function of the number of independent samples, N. For the present digital 
computation. 


N = 65,536 
e [ AP(RMS) ] = 

Thus, for both analog and digital computations, the normalized standard error in the 
overall RMS fluctuating pressures is less than 1 percent and may be disregarded as a 
significant uncertainty. 

To assess the technique inaccuracies in the overall RMS fluctuating pressure level, it 
suffices to use a comparative technique whereby the analog computations are com- 
pared with the digital computations. The basis for this analysis is that technique 
errors and uncertainties are primarily associated with the analog instrumentation and 
the personnel performing the data reduction; whereas, the digital computations are 
relatively free of technique errors and uncertainties. Noting that both analog and 
digital computations are free of significant statistical uncertainties, it suffices to 
evaluate the technique inaccuracies from a comparison of the analog and digital 
computations. Data, representing 286 calculations for the Mach number range from 
0.60 to 1 .60, are presented in Figure 20. It is evident from the results presented in 
Figure 20 that the analog results are, in general, greater than the digitally computed 
results. Therefore, it appears that the analog results have been biased to the high 
sick;, although this error is very slight - on the order of 4 percent or 0.4 dB. This 
error is well within the accuracy of the analog instrumentation employed in the data 
reduction process. The scatter about the mean, which is representative of the tech- 
nique uncertainty in the data, is bounded by ± 8 percent with 95 percent confidence. 
That is, 95 percent of the data falls within ± 8 percent of the mean. Thus, assuming 
a Gaussian distribution of the scatter about the mean, the normalized standard error - 
which is based on a 68 percent confidence level - is one half of the 95 percent con- 
fidence interval, or e = 0.04. In summary, the technique error is given by: 

e [AP(RMS)] = 0.04 [AP(RMS)] (5.24) 




2N = 0 - 0026 


(5.23) 
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and the technique uncertainty is given by 


± e [AP(RMS)] = ±0.04 f AP(RMS) ] (5.25) 

In coefficient form (neglecting the uncertainty in q ), Equations 5.24 and 5.25 
become 00 

e [AC p (RMS)] t = 0.04 [AC p (RMS)] (5.26) 

± e [AC (RMS)] = ±0.04 [AC (RMS)] (5.27) 

r » r 


Power Spectral Density 

Uncertainties in the power spectra are dominated by the statistical uncertainties since 
these data were reduced to descriptive terms by digital computation. It is felt that 
technique errors and uncertainties are very small for the power spectra since calibra- 
tions recorded on tape were used for the specification of the sensitivity for each 
microphone without any intermediate conversion or interpretation of the data by 
personnel . Certainly, the technique errors are negligibly small compared to the 
statistical uncertainties. The normalized standard error for conventional digital 
computation of the power spectral density is given by 

e t<f(f)I s = Vm/N (5.28) 

where m is the maximum number of correlation lag values of the computation and N is 
the size of the digital sample. During the present computations, a rather complicated 
digital filtering technique was employed whereby a series of power spectral density 
points were calculated with each having 16 degrees of freedom. The method of cal- 
culating the power spectral density points involved dividing the frequency range into 
a series of frequency intervals such that the power spectral density for each interval 
was computed with 16 degrees of freedom. This technique is discussed further in 
Section 4.4.2. The power spectral density for each 1/3 octave frequency was then 
computed by averaging the power spectral density for the number of discrete points (k) 
within the 1/3 octave band. The resulting number of degrees of freedom for each 
1/3 octave band thus becomes 16 k. To increase the statistical accuracy at the lower 
frequencies, the 1/3 octave data were further averaged to give octave band results 
for the octave band center frequencies of 31 .5 Hz, 63 Hz and 125 Hz. The resulting 
variation of the statistical number of degrees of freedom, n, with frequency is pre- 
sented in Figure 21 . The normalized standard error for the power spectral density may 
be determined directly from the number of degrees of freedom as follows: 


(5.29) 


* (f(f)L 

o 



The resulting variation of e [ ] with frequency is presented also in Figure 21 . 

s 

5.3.4 Summary of Errors and Uncertainties 

The following is a summary of the errors and uncertainties associated with the test 
results presented herein. The errors and uncertainties are presented in equation form 
with the required terms referenced to figures presented in the report where necessary. 

Overall Level 

Error 

e [ AC p (RMS) ] = {e [ A C p (RMS) ] | + e f AC p (RMS) ] R + e [ AC p (RMS) ] 2 }* 

1 

= {e [ ACp(RMS) I 2 e+r + e [AC p (RMS)J 2 } 2 

where 

e (AC p (RMS)J 2 +R = [ AC p (RMS) ] 

e [ AC (RMS) J 2 = 0.0016 [AC (RMS)] 2 

Pt P 

e [<|>(f)] £ * s given in Figure 16, and 

R(f) is given in Figure 19. 

Uncertainty 

± e f AC (RMS) J = db e [ACJRMS)]. 

P P I 

= ±0.04 f AC (RMS)] 

r 
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Power Spectral Density 


Error 

e (cj>(f)] = e [ <|>(f) ] E + e [$(f)] R 


6 [ < j >(0 ] 

|*H0 - e [q>(f)] E j” 

W) 


= *(f) - 


Uncertainty 

± € [<}>(f)] = ± € [ <t>(0 ] $ 

where e [ <|>(f ) ] is given in Figure 21 . 
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6.0 


DISCUSSION OF TEST RESULTS 


6.1 Introduction 

Because of the large number of data measurements obtained during the test, it is 
necessary to be selective in the presentation and discussion of the test results. Noting 
that the primary objective of the present report is to define the basic features of the 
flow field induced by three-dimensional protuberances, the measurements obtained 
along the longitudinal centerline of the model were selected as being primary data. 
These data represent the plane of symmetry of the flow around the protuberances and, 
thus, show the basic features of the flow in both the upstream separated region and 
the downstream wake region. The flow around the protuberances, which is charac- 
terized by the transition from separated to wake flow is discussed only briefly. A 
comprehensive presentation of the data for the flow around the protuberances will be 
made in a later report. 

The test results are discussed under the following section titles: 

• Clean Model Flow Characteristics 

• Protuberance Induced Flow Characteristics 

— Protuberance Induced Static Pressure Field 
— Protuberance Induced Fluctuating Pressure Field 

• Analysis of the Protuberance Induced Flow Field 

The clean model flow characteristics were determined from the test panel measurements 
for the protuberance in the fully retracted position. Results for the clean model indi- 
cate the quality of the flow over the test panel and the suitability of the mode! for the 
investigation of protuberance flow fields. The static pressure field induced by the 
protuberances is exhibited by surface pressure measurements taken over the test panel 
and the wall of the protuberance. Similarly, the fluctuating pressure field induced by 
the protuberance is exhibited by surface microphone measurements. Finally, the static 
and fluctuating pressure measurements are combined with the oil-flow visualization 
results to determine the salient features of the protuberance induced flow field. 

6.2 Clean Mode! Flow Characteristics 

6,2.1 Surface Pressure Distributions 

Static pressure data as obtained from the distribution of orifices along the longitudinal 
and transverse centerlines of the test panel are presented in terms of the pressure 
coefficient, Cp, in Figure 22. The variations of C over the surface of the test panel 
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are significantly greater than the uncertainty associated with the measuring instru- 
mentation and, also, these data show some consistent trends. Thus, it is concluded 
that the variations are indicative of the aerodynamic characteristics of the model, 
or the model/tunnel combination. The longitudinal variations in C , Figure 22a, 

exhibit the irregularities of the clean model flow field, and when compared with the 
transverse variations. Figure 22b, it is evident that the flow is less uniform in the 
longitudinal direction than in the transverse direction. It is interesting to note the 
agreement between the longitudinal and transverse pressure distributions as indicated 
by the solid symbols on each figure. The agreement is within the accuracy of the 
instrumentation employed. This further substantiates the conclusion that static pres- 
sure variations result from aerodynamic effects associated with either the model or the 
model/tunnel combination. Since the variations appear to be greater at > 1 .0 

than at subsonic Mach numbers, it is felt that extraneous Mach waves in the tunnel 
contribute significantly to the static pressure variations. Also, the boat tail of the 
model splitter-plate appears to have affected the longitudinal centerline pressures at 
Mqq = 1 .20. In general, the static pressure variations about the free-stream pressure 

are within ±0.1 of the free-stream dynamic pressure. Normally, variations should be 
contained within the uncertainty limits of the instrumentation. Thus, it is felt that 
the quality of the flow over the test article was not as good as desired; however, the 
flow is acceptable for the study of protuberance flow fields since the clean model 
variations in pressure are relatively small compared with the pressure variations 
induced by the protuberances. 


At this point it should be noted that the static pressure results for the protuberance 
induced environments were adjusted for the clean model results by subtracting the 
static pressures for the protuberance-free case from the static pressure for the 
perturbed case. This adjustment contributes to the uncertainty of the final results; 
however, a significant improvement in the static pressure distributions was obtained 
for the perturbed flow. Thus, it is felt that this adjustment was justified. 


6.2.2 Boundary Layer Characteristics 

Further information of the charac ('eristics of the flow over the clean model may be 
gained from an examination of the boundary layer characteristics. Boundary-layer 
velocity profiles were determined from pitot pressure measurements obtained with the 
forward and aft rakes with the protuberance in the fully retracted position. The 
velocity profiles, at the locations of the forward and aft rake, are presented in 
Figure 23. Since the rakes were positioned near the leading and trailing edges of the 
test panel at approximately equidistance from the location of the protuberance, the 
boundary layer characteristics at the protuberance location may be taken as the mean 
of the data for the forward and aft rakes. 

It is evident from the test results that an extremely thick boundary layer existed over 
the region of the test panel for M^ = 0.60 and 0.80. A similar environment was 

observed in the small scale tests reported in Reference 18. This environment is 
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attributed to a partially blocked secondary flow under the splitter plate which induces 
local separation at the leading edge of the model . The separation and subsequent 
reattachment of the flow at the leading edge results in a thickened boundary layer 
over the test panel. Also, it will be noted that the velocity ratios do not approach 
1 .0 at the edge of the boundary layer at every Mach number. This anomaly resulted 
from the deviation of the local velocity from free-stream conditions as indicated by 
the variations in static pressure presented in the preceding section. 

The variations of the boundary layer displacement thickness, 5*, with free-stream 
Mach number are presented in Figure 24. The displacement thicknesses for the rake 
data were determined from numerical integrations of the velocity profiles. For free- 
stream Mach numbers of 0.90 and above, the displacement thickness varied irregularly 
with Mach number about an approximately constant value. For Mqq = 0.60 and 0.80, 
the thick boundary layer resulting from leading edge separation is clearly evident in 
the large values of S*. The displacement thicknesses shown for Mqq = 0.60 and 0.80 

were obtained by numerical integrations to the point where U n is approximately 0.95 

Uqq, thus, it is anticipated that 5* may be slightly greater than the values shown. 

The experimental results are compared with an empirical prediction using Bies' formula 
(see Reference 22). The empirical prediction represents the estimated displacement 
thickness at the location of the protuberance which would correspond to the mean 
between the forward and aft rake results. The empirical results agree closely with the 
results presented for the forward rake, indicating that the boundary layer at the pro- 
tuberance location was somewhat thicker than predicted. Surface excrescences in the 
boundary layer upstream of the test panel apparently caused a slight artificial thickening 
of the boundary layer. 

The variations of boundary layer momentum thickness, 0, with free-stream Mach 
number are shown in Figure 25. For the Mach number range of the present test, 
momentum thickness is essentially invariant with Mach number except for the anomaly 
due to the separation induced thickening of the boundary layer at M^ = 0.60 and 

0.80. Again, agreement is shown between the experimental measurements for the 
forward rake and the empirical results. 

The effects of free-stream unit Reynolds number on the boundary layer characteristics 
in the vicinity of the test panel (Figures 24 and 25) were negligibly small . This is as 
anticipated since, for a turbulent boundary layer on a flat plate, the boundary layer 
thickness is inversely proportional to the Reynolds number to the one-fifth power. 

A summary of important boundary layer and free-stream flow conditions which are 
characteristic of the undisturbed flow over the model in the vicinity of the protu- 
berance for a unit Reynolds number of 3 x 10 6 are presented in the following 
table . 


UNDISTURBED FLOW PARAMETERS OVER THE TEST PANEL 

Re/ ft = 3x 10 6 /ft 


M„ 

00 

Uoo 

ft/sec 

psf 

mm 

mm 

0 

in . 

0.60 

676 

388 

«8 

» 0.80 

« 0.50 

0.80 

880 

479 

« 8 

« 0.80 

» 0.50 

1 .00 

1067 

559 

2.5 

0.36 

0.21 

1 .20 

1264 

611 

2.5 

0.37 

0.19 

1 .40 

1391 

644 

2.5 

0.36 

0.19 

1 .60 

1528 

662 

2.5 

0.31 

0.20 


6.3 Protuberance Induced Flow Characteristics 


Before discussing the static and fluctuating pressure measurements, some general 
comments in regard to the protuberance flow field should be made. A photograph 
of a typical oil flow pattern around the protuberance is presented in Figure 26. Any 
protuberance which extends into or through an attached boundary layer will (I) tend 
to induce separated flow upstream of the protuberance and (2) generate a wake down- 
stream of the protuberance. For a three-dimensional protuberance the upstream 
separated flow will have a curved separation line or front as clearly shown in Figure 26. 
The wake will be characterized by a short converging region which is terminated by a 
narrow neck followed by a diverging region. The converging region will be denoted 
as the near wake and the diverging region as the far wake with the neck region 
separating the two. The neck and far wake are clearly shown in Figure 26; however, 
the near wake is partially obscured by the protuberance due to the slight downstream 
angle of the camera view. There are several other interesting features of the flow 
field shown in Figure 26; however, these will be discussed in detail in later sections. 

At this point, it suffices to recognize only the presence of the upstream separated flow 
and the downstream wake regions, and that the region of perturbed flow may extend 
for several diameters around the protuberance. 

The length of the upstream separated flow region at the longitudinal centerline was 
evaluated from oil-flow studies of the 4-inch diameter protuberance for the Mach 
number range from 0.60 to 1.40 and h/ D range from 0 to 2.0. The separation point 
is identified in Figure 26. The resulting variations of normalized separation length, 
i s /D, with h/ D for the various Mach numbers of the test are presented in Figure 27. 

The curve presented for Mqq = 1 .60 was determined by extrapolation of the lower Mach 
number data to Mqq = 1 .60. An increase in either protuberance height or Mach number 

resulted in an increase in separation length. The trend in the data indicates that the 
separation length will approach asymptotic values for large protuberance heights. 
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In the wake, it is obvious that the flow will be perturbed for many diameters down- 
stream of the protuberance. The aft rake, which was present throughout the test, 
revealed noticeable distortion in the velocity profiles due to the wake for the various 
protuberances. For the 2-inch diameter protuberance, the rake was located at 18.28 
diameters aft of the protuberance and the effect of the wake was barely detectable. 
Thus, it is felt that the test panel surface measurements for the three protuberance 
configurations adequately cover the region of perturbed flow as discussed in following 
sections. Generally, the neck region which separates the near- and far-wake regions 
appeared to occur at approximately one diameter aft of the protuberance (r/ D « 1 .5) 
for Mqq > 1 .20. For lower Mach numbers, the wake was ill defined by the oil flow 

pictures; however, the fluctuating pressure measurements indicate that the neck region 
moved aft with decreasing Mach number. 

6.3,1 Protuberance Induced Static Pressure Field 

The static pressure field induced by the protuberances is represented by the distribution 
of the pressure coefficients along the test panel centerline, forward and aft of the pro- 
tuberance, Figures 28-30. Test results are presented, separately, for the 2-, 4-, and 
8-inch diameter protuberances for various heights ranging up to two diameters. Tlv: 
2-inch diameter protuberance results. Figure 28, show the extent of the perturbed fiow 
field upstream and downstream of the protuberance; however, little detail of the pres- 
sure field is given for the region in close proximity to the protuberances. The 4-inch 
diameter protuberance results. Figure 29, give a cursory definition of the perturbed 
flow field for the region in close proximity to the protuberance as well as the outer 
extremities of the protuberance environment. Finally, the 8-inch diameter protuber- 
ance results. Figure 30, show, in detail, the perturbed flow field for the region in 
close proximity to the protuberance. 

An examination of the results presented in Figures 28-30 reveals several distinct 
characteristics of the pressure field which are typical of cylindrical protuberances. 

The static pressure field associated with the upstream separated flow is characterized 
by positive pressure coefficients; whereas, the static pressure field associated with 
the wake is characterized by 1) negative pressure coefficients in the near wake and 
neck regions, and 2) approximately zero pressure coefficients in the far wake. Clearly, 
the regions most affected by the protuberance lie in close proximity to the protuber- 
ance. It is important to note that the region of perturbed flow upstream of the 
protuberance is not bounded by the region of separated flow, but may extend beyond 
this point. A comparison of the separation points (Figure 27) with the static pressure 
distributions will show this characteristic. At subsonic and slightly supersonic Mach 
numbers, the perturbed flow field extends significantly beyond the separation point; 
whereas, at M^ > 1 .20, the shock strength is sufficiently large to cause separation 

near the foot of the shock wave. For the latter Mach number range, the region of 
upstream perturbed flow (which is bounded by the shock wave) and the region of 
separated flow occur over approximately the same area with the separated flow region 
being contained slightly within the perturbed flow region. Consistent with this finding 
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are the trends in the pressure gradients associated with the upstream perturbed flow. 

At subsonic Mach numbers, a gradual increase in siatic pressure in the flow direction 
is observed. The effect of the positive pressure gradient is to retard the boundary 
layer flow near the wall - eventually causing separation. The critical pressure at 
separation is discussed below. At supersonic Mach numbers, a rapid increase in static 
pressure is observed across the shock wave; thus inducing separation immediately aft 
of the mean location of the shock wave. The variations of the static pressure ratio at 
separation with Mach number are presented in Figure 31 . The results for the different 
size protuberances were obtained based on the separation lengths as defined from oil 
flow patterns for the 4-inch diameter protuberance. It should be noted that separation 
length probably varies with protuberance diameter rather than separation pressure such 
that all three configurations would have the same separation pressure ratio as given by 
the 4-inch diameter protuberance data. 

Downstream of the separation point, a region of relatively constant, plateau, pressure 
occurs. The plateau pressure region is particularly noticable for M^ > 1 .20. For 

two-dimensional separation, the plateau pressure extends to near the step face; how- 
ever, for the three-dimensional case, the plateau region is terminated significantly 
upstream of the protuberance. For the latter case, the plateau region is terminated 
by a sudden reduction in static pressure which is followed by a sharp increase in 
pressure near the protuberance. The. low pressure region imbedded within the 
separated flow field is unique to the three-dimensional protuberance case and pro- 
vides much insight into the understanding of the structure of the three-dimensional 
separated flow field. The separated flow field is discussed at length in Section 6.4. 

The wake region exhibits noticable variations in static pressure for approximately 
four diameters aft of the protuberance. It is interesting to note that distortion in the 
velocity profiles were observed for approximately 18 diameters aft of the protuberances 
as previously mentioned; however, the retarded flow in the wake apparently has little 
effect on the surface static pressures in the far wake region. Little insight into the 
structure of the wake can be gained from the static pressure measurements and it is 
felt that the oil-flow photographs and fluctuating pressure measurements are more 
definitive of the wake structure. 

A cursory examination of the effects of Reynolds number and protuberance size on the 
perturbed flow field was made utilizing the static pressure measurements. The effects 
of Reynolds number were evaluated using the 8-inch diameter protuberance results 
for unit Reynolds number conditions of 1 .5, 3.0 and 4.5 million per foot. Corre- 
sponding Reynolds numbers, based on protuberance diameter, are 1 .0, 2.0, and 3.0 
million. Representative axial distributions of static pressure coefficient are pre- 
sented for h/ D = 2.0 and Mach numbers from 0.60 to 1 .60 in Figure 32. The effect 
of Reynolds number variations appear to be small, although there are some noticable 
deviations in the data. In general, variations in Reynolds number primarily affected 
the static pressure distributions in the wake region at subsonic Mach numbers (although 
a consistent trend was not evident); whereas, at M^ > 1 .20, the effects appeared to 
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be more predominant upstream of the protuberance. Increasing the Reynolds number 
at supersonic Mach numbers appeared to cause a slight enlargement of the upstream 
separated flow region. The effects of protuberance size variations were determined 
by comparing axial distributions of static pressure coefficient for the 2- and 4-inch 
diameter protuberances with Reynolds number, based on protuberance diameter, held 
constant at 0.5 million (Figure 33). Further, a similar comparison was made between 
the 4- and 8-inch diameter protuberances for a Reynolds number, based on protu- 
berance diameter, of 1 .0 million (Figure 34) . These data represent the h/D = 2.0 
condition for Mach numbers ranging from 0.60 to 1 .60. Small variations in static 
pressure coefficient were observed in both the upstream and wake regions due to pro- 
tuberance size variations with the most significant effects occurring in the wake 
region at subsonic Mach numbers and in the upstream separated flow region at super- 
sonic Mach numbers. The most noticeable effect of increasing the protuberance size 
was to cause a slight enlargement of the normalized upstream separated flow region 
at supersonic Mach numbers (particularly for Rep = 0.5 x 10 6 ). 

The pressure distributions over the wall of the 8-inch diameter protuberance are 
exhibited by variations of Cp with vertical location, y /D for various polar angles 

around the protuberance. Test results are presented for a range of protuberance 
heights and Mach numbers in Figures 35-40. A high density of static pressure orifices 
were present along the vertical rays and thus provided good definition of the pressure 
variations associated with the various flow fields to which the protuberance was 
exposed. Flow fields impinging on the upstream half of the protuberance (0 + < 90 
degrees) which appear to be evident in the static pressure measurements are: 1) the 
separated shear layer (base region) near the base of the protuberance, 2) the two- 
dimensional cross flow (central region) near the middle of the protuberance, and 
3) the end effect (tip region) near the top of the protuberance. The vertical extent 
of these regions naturally vary with h/D. For h/ D > 1 .5, all three regions are 
clearly defined in the vertical pressure distributions. The tip region is characterized 
by a combination of flow around the protuberance and over the top of the protuberance. 
The resulting three-dimensional relief afforded to the flow causes a reduction in pres- 
sure in comparison to the pressures for two-dimensional flow as encountered by an 
infinite cylinder in cross flow. At supersonic Mach numbers, the bow shock generated 
by the cylinder is curved in the vicinity of the tip region as discussed in References 17 
and 18. The tip region extends downward from the top of the protuberance for approxi- 
mately 0.4 diameters for h/D,,* 1.5. The central region is characterized by nearly 
parallel or two-dimensional flo>v around the protuberance. In the central region, the 
wall pressure is approximately constant along vertical rays and the flow exhibits 
characteritistics similar to that of an infinite cylinder in cross flow. The central 
region is located between the rip region and the base region and its vertical extent 
varies noticeably with h/D. The base region is characterized by the interaction of the 
separated flow, which is induced upstream of the protuberance, with the protuberance. 
The presence of reverse flow in the separated flow region is evidenced by the low pres- 
sures measured just above the base of the protuberance. For Mqq > 1 .0, the growth of 
the separated flow region is exhibited by the upward movement of the pressure defect 
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in the vertical pressure distributions. For Mq-, < 1 .0, the pressure defect associated 

with the base flow is poorly defined due to the small extent of the separated flow 

region. Further, the region of constant pressure normally associated with the central 

region is not exhibited in the results for Mqq <1.0 due to the extremely thick inflowing 

boundary layer for M^ = 0.60 and 0.80. The boundary layer velocity profile shape in 

the region of y /D < 1 .0 results in a loss of total pressure and, thus, a decrease in the 
P 

surface pressures measured over the front half of the cylinder. 


Over the aft half of the cylinder (0 + > 90 degrees) small variations in pressure over 
the wall of the cylinder were observed for M 00 > 0.60. This is explained by the 
separation of the cylinder boundary layer near the 0 = 90-degree ray. This charac- 
teristic is clearly evident in the test results for M^ >1.0. For a Mach number of 

0.60, large variations in the protuberance wall pressures were observed over the wall 
of the cylinder at the 0 + = 90-degree ray. Apparently the flow exhibits irregular 
separation along the 90-degree ray due to the variation in the flow upstream of the 
protuberance in the vertical direction. 


6.3.2 Protuberance Induced Fluctuating Pressure Field 
Overall Levels 

Of particular interest during the present study was the fluctuating pressure field 
induced by the cylindrical protuberances. The overall fluctuating pressure levels 
associated with the perturbed flow field are presented as variations of the RMS 
coefficient of fluctuating pressure along the test panel longitudinal centerline in 
Figures 41-43. Test results are presented, separately, for the 2-, 4-, and 8-inch 
diameter protuberances for ranges of protuberance height and Mach number which are 
consistent with the presentation :>f static pressure results discussed in the preceding 
section. Here again, the results for the 2-, 4-, and 8-inch diameter protuberances 
represent varying degrees of definition of the perturbed flow field ranging from a 
general overview as exhibited by the 2-inch diameter protuberance results to a 
detailed, near field, study as exhibited by the 8-inch diameter protuberance results. 

Certain characteristics of the fluctuating pressure field are immediately obvious from 
an examination of these results. The fluctuating pressure field induced by the protu- 
berances corresponds closely with the region of perturbed flow as exhibited by the 
static pressure results. Two distinct regions of perturbed flow are exhibited by the 
fluctuating pressure results - one corresponding to the upstream separated flow region 
and the other corresponding to the downstream wake region. It is evident from the 
fluctuating pressure measurements that significant fluctuating pressures occurred in 
both the upstream and downstream environments. In general, the most severe fluctuat- 
ing pressures were encountered in the wake region for M 0 0 < 1 .0 and large h/ D; 
whereas, for other test conditions, the separated flow upstream of the protuberance 
coused the most severe fluctuating pressures. It is convenient to discuss the separated 
and wake flow regions separately. 
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In the separated flow region, the fluctuating pressure field is characterized by a 
gradual increase in overall level as the field is observed in the direction of the 
free-stream. Near the mid-point in the separated flow region, the fluctuating 
pressure reaches a maximum value which ranged from 12 to 17 percent of the free- 
stream dynamic pressure. The peak is followed, in general, by a decrease in 
fluctuating pressure level . This decrease appears to become less drastic with 
increasing Mach number such that, at = 1 .60, a peak plateau level of fluctuating 

pressure occurs near the protuberance. For M^ =1.40 and 1 .60, a small peak in 

fluctuating pressure occurs upstream of the separation point and corresponds to the 
mean location of the shock wave. Levels of fluctuating pressure for the shock wave 
are approximately the same as observed for two-dimensional protuberances (forward 
facing steps); however, the levels within the separated flow region are considerably 
greater than for the two-dimensional case. Immediately, it becomes obvious that 
the unsteady mechanisms within the three-dimensional separated region are significantly 
different than the two-dimensional case and, further, the three-dimensional flow 
generates more severe fluctuating pressures. 

Increasing either protuberance height or free-stream Mach number results in increasing 
the area of the unsteady perturbed flow field upstream of the protuberance. For the 
. range of test variables of the present study, increasing the protuberance height did not 
alter the peak level of the fluctuating pressures in the upstream separated flow region 
with any consistent trend and, in general, there appears to be only a slight increase 
in the peak level of the fluctuating pressure with increasing Mach number. However, 
the most significant and important trend is the increase in area which is subjected to 
the high levels of fluctuating pressure within the separated flow region for increasing 
protuberance height and free-stream Mach number. 

The fluctuating pressure field associated with the wake of the protuberance extended 
over a relatively large axial distance for all protuberance heights and Mach numbers 
as is shown in Figures 41-43. Very large fluctuating pressures were encountered at 
subsonic Mach numbers, and in certain cases, the levels exceeded those encountered 
upstream of the protuberance. Of particular interest are the results for M^ = 0.80 

and h/D = 1 .5 and 2.0, Figure 43b. For these conditions, the peak levels in the 
wake ranged from 18 to 22 percent of the free-stream dynamic pressure . For the unit 
Reynolds number of 3 x 10 6 , the corresponding overall sound pressure level ranges 
from 166 to 168 dB. These values were the largest encountered in the perturbed flow 
field along the longitudinal centerline of the test panel for the range of test variables 
investigated. As previously mentioned, it was difficult to distinguish the neck region 
of the wake at Mgg < 1 .20; however, there are indications that the neck occurs further 

downstream than the location of the peak levels of fluctuating pressure for Mq^ = 0.80. 

(Some evidence of this is given by the cross-spectra results as discussed in later para- 
graphs.) For certain conditions (for example, Mqq = 1 .00, h/ D = 2.0) the peak levels 

in the fluctuating pressures associated with the wake occurred 3 or 4 diameters aft of 
the protuberance. It is difficult to define consistent trends in the wake results, and 
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perhaps this confusion is a result of the complex flow structure that exists in the wake 
and the way in which it changes with variations in both Mach number and protuberance 
* height. Additional discussion on this point is presented in the section on the power 
spectra . 

The general trend in the axial distributions of the fluctuating pressures in the wake 
consists of a moderate increase to peak levels followed by a gradual decrease to the 
unperturbed flow level as the flow field is examined in the direction of the free- 
stream. Both the location and the width of the peak varies with free-stream Mach 
number and protuberance height. Further, for M^ > 1 .0, there appeared to be a 

general reduction in the overall fluctuating pressures in the wake with increasing 
Mach numbers. Also a gradual aft movement and decrease in level of the peak 
fluctuating pressures occurred with increasing protuberance height. 

A comparison of the axial distributions of overall fluctuating pressure level for the 
2-, 4- , and 8-inch diameter protuberances is shown in Figure 44. All of the pro- 
tuberances exhibited similar trends in the axial distributions; however, some devia- 
tions in the levels were noted. These deviations are relatively small and may result 
from either Reynolds number (based on protuberance diameter) or protuberance size 
effects . 

The fluctuating pressure field over the surface of the protuberance is exhibited by 
measurements taken over the wall of the 8-inch diameter protuberance. Figures 45- 
50. These data are presented as variations of RMS fluctuating pressure coefficient 
along circumferential rays with 0" = 0 degree corresponding to the upstream center- 
line and 0“ = 180 degrees corresponding to the downstream centerline. It will be 
noted that a range of circumferential rays are presented for various heights above 
the test panel surface. For h/ D = 1 .0, the surface measurements were contained 
within the perturbed flow field of the test panel (base region); whereas, for h/ D = 

2.0, the surface measurements were above the perturbed flow field of the test panel 
and should closely approximate the characteristics of an infinite cylinder in cross 
flow (central region). In general, flow around the cylinder is characterized by 
attached flow over approximately the forward half of the cylinder and separated flow 
over the aft half. The separation of the flow from the wall of the cylinder occurred 
near the 90-degree polar location. For M^ = 0.60, h/D = 2.0, extremely large 

fluctuating pressures (in excess of 170 dB) resulted from the unsteadiness of the separa- 
tion line (see Figure 45a). These data were rechecked following the initial discovery 
of the large levels to verify their accuracy . Further, it should be noted that these 
results are consistent with the location of the separation line as observed in the static 
pressure results. Some evidence of flow separation near 0" = 90 degrees is given by 
the peak in the circumferential distributions at this location for all M ^ < 1 .00; 

whereas, for M^ > 1 .00, it is felt that separation occurs slightly downstream of the 

90-degree location. Due to the present model design, a continuous survey of the 
polar distribution could not be made. For h/D = 2.0, little variation with vertical 
location was observed in the fluctuating pressure levels except for the = 0.60 
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condition. Further, the wake region (9” > 90 degrees) appeared to exhibit larger 
fluctuating pressure levels than did the upstream region at subsonic Mach numbers. 

For h/D = 1 .0, the effect of the perturbed flow field around the base of the protu- 
berance is evidenced by the large fluctuating pressures in this region. For Mqq > 1 .20, 
the fluctuating pressure levels over the front half of the cylinder appeared to decrease 
with an increase in polar angle. It is interesting to note the drastic reduction in 
fluctuating pressures for the rear half of the protuberance at M^ > 1 .20. These data 

appear to be consistent with an overall change in the wake structure at supersonic 
Mach numbers as observed in the oil flow studies and static and fluctuating pressure 
measurements taken over the test panel. Further, the test panel data for microphones 
near the protuberance appear to agree reasonably well with the protuberance data at 
0” = 0 degree and 180 degrees as shown in Figures 45-50. 

Power Spectra 

Power spectra represent the distribution of the mean square energy of the fluctuating 
pressure with frequency. In general, power spectra of most fluctuating pressure 
environments are found to scale on a Strouhal number basis; that is, the frequency 
is non-dimensionalized by multiplying a typical length and dividing by a typical 
velocity. The advantages of using normalized spectra are obvious since it enables 
similar, homogeneous, flow fields to be represented by a single spectrum regardless 
of the scale of the flow field or the free-stream velocity. Numerous studies have 
been made to determine the proper parameters to be used to non-dimensional ize the 
spectra for various aero-acoustic environments. Unfortunately, the choice of 
parameters which best collapses the data appear to be dependent on the nature of 
the fluctuating pressure environment. In general, free-stream velocity is used as 
the normalizing velocity parameter, although a typical eddy convection velocity 
(itself a function of frequency) has been used occasionally. The local convection 
velocity appears to correspond more closely with the physical situation for fluctuating 
pressures due to turbulent eddies. For simplicity, free-stream velocity has been 
adopted for the power spectra presented herein since the three-dimensional protuber- 
ance flow field is a non-homogeneous environment with widely varying convection 
velocity characteristics. 

Selection of a typical length is more difficult. For both undisturbed boundary layer 
flow and two-dimensional separated flow, typical dimensions of the undisturbed 
boundary layer such as boundary layer thickness, displacement thickness, and 
momentum thickness, have all been used by various authors. For transonic flow, the 
typical boundary layer lengths vary only slightly, but no final conclusion can be 
drawn on the relative merits of the collapse for any particular length. The most 
generally used typical length is the boundary layer displacement length, 5*, and 
this parameter is used herein. An important typical length parameter for the pro- 
tuberance flow field is the protuberance diameter, particularly for the region of the 
wake. It is evident from the test results presented herein that the scale of the perturbed 
flow field varies in proportion to the protuberance diameter. Thus it must be concluded 
that the power spectra of the fluctuating pressures would depend to some degree on the 
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size of the protuberance. However, since test results are presented for the 8-inch 
diameter protuberance only, the use of protuberance diameter as a typical length 
parameter will be withheld until the power spectra for the other size protuberances 
are reduced. These data will be forthcoming in a later report. 

Power spectra representing the frequency characteristics of the protuberance induced 
fluctuating pressures are presented for both upstream and downstream flow fields in 
Figures 51-62. It will be noted that test results are presented for the Mach number 
range from 0.60 to 1 .60 and protuberance heights of 1 .0 and 2.0 diameters. These 
data were obtained for the 8-inch diameter protuberance configuration and, thus, pro- 
vide a detailed study of the protuberance induced flow field for the region in close 
proximity to the protuberance. It is impractical to discuss all of the frequency charac- 
teristics of the flow as exhibited by the power spectra; however, some general features 
should be pointed out. The data for the upstream and downstream flow fields will be 
discussed separately. 

For the upstream flow field, test results are presented only for the microphone locations 
which appeared to have been within the perturbed region (Figures 51-56). All of the 
spectra appear to exhibit the same general characteristics when presented in terms of 
normalized power spectral density. For the present discussion, it suffices to note that 
a high concentration of energy was present at low frequencies with a roll-off at the 
high frequencies. For the region of the most intense fluctuating pressures upstream 
of the protuberance at < 1 .0, the roll-off appeared to be approximately two 

orders of magnitude per decade. This is equivalent to 20 dB per decade or 6 dB per 
octave band. For < 1 .0, there appears to be relatively little variation in the 

slidpe of the spectra within the separated flow region, although there is some change 
from the perturbed flow region upstream of separation to the region within the 
separated flow field. For the perturbed flow field upstream of the separation point, 
the spectra appear to roll-off less rapidly at the high frequencies. For > 1 .0, 

the separated flow field is sufficiently large in comparison with the microphone spacing 
to enable a trend in spatial variation of the power spectra to be established. Of par- 
ticular interest is the plateau in the power spectra in the vicinity of r/D = - 1 .00 to 

- 1 .25 for h/D = 1 .0, Mqo = 1 .40 and 1 .60 and in the vicinity of r/ D = - 1 .125 to 

- 1 .500 for h/D = 2.0, = 1 .40 and 1 .60. This region corresponds to the inner 

separated flow field in the vicinity of the peak levels of fluctuating pressure (see 
Section 6.4). The plateau indicates un increase in fluctuating pressure energy at 
the mid and high frequencies when compared to the spectra for other regions of the 
separated flow field. 

Another important characteristic of the perturbed flow field is exhibited by the power 
spectra associated with the shock wave at = 1 .40 and 1 .60, h/D = 1 .0. At the 

mean location of the shock wave (r/D = - 1 .25), the flow is disturbed as a result of 
the shock wave impinging on, and interacting with the attached boundary layer. The 
power spectra for these data (Figures 55a and 56a) exhibited a sharp reduction in 
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energy with increasing frequency starting at relatively low frequencies. The 
transition in spectra shape is evident also as the fluctuating pressure field is 
examined moving aft of the shock wave . The transition is characterized by a 
reduction in low frequency energy and an increase in high frequency energy. 

Clearly, the power spectra show a large spatial variation in the fluctuating pressure 
characteristics as was observed in the overall fluctuating pressure levels. The upstream 
perturbed flow field is obviously non-homogeneous in space in both overall level and 
frequency characteristics. However, it is interesting to note that variations in protu- 
berance height and free-stream Mach number did not appear to have any gross effect 
on the general shape of the power spectra for a given region in the perturbed flow 
field. Although a direct comparison is not presented, a cursory examination of the 
data suggests that spectra for the same values of r/£ g would agree at least in trend 

for the various h/D and conditions. This characteristic will be evaluated in 

detail in a subsequent report and, if correct, it will be a strong argument for the 
similarity of the upstream separated flow fields for various Mach numbers and 
protuberance heights. 

The power spectra for the wake. Figures 57-62 exhibit the characteristics of a non- 
homogeneous flow field which is gradually changing in space. It is difficult to 
determine precise trends in the data with variations in either Mach number, protu- 
berance height, or axial location. The general trend is that the spectra appear to 
roll-off less rapidly at the higher Mach numbers. Also, the spectra for supersonic 
Mach numbers have appreciably less low frequency energy than the subsonic data. 

One interesting feature of the wake spectra is the peak which occurs at certain 
conditions. For example, the spectra for h/D = 1 .0, bA^ = 1 .20 (Figure 60a) show 

a peak at f S*/ =0.013 for all axial locations. The corresponding Strouhal 

number, based on protuberance diameter, is 0.28. The peak generally occurs at the 
same Strouhal number for other conditions; however, the peak does not occur at all 
conditions. Further, the value of 0.28 is larger than the Strouhal number for vortex 
shedding from an infinite cylinder in cross flow which has a value of 0.21 for the 
Reynolds number range of the present test. However, it is probable that the peaks 
observed in the present spectra are due to vortex shedding from the cylinder at a 
Strouhal frequency of 0.28. 

Cross-Spectra 

Considerable insight into the unsteady mechanisms at work in a fluctuating pressure 
field can be gained by studying the relationship between data observed at two dif- 
ferent points in the field. For a number of years, investigators involved in turbulence 
and unsteady-fluid dynamic research have utilized cross-correlation techniques to 
describe, on the average, the mutual statistical dependence of the signals at two 
different points. From such studies, the correlation patterns, convection velocities, 
and other features of the flow have been evaluated. One useful technique which 
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has been adopted in recent years is the cross* 5 * spectra I density function. The cross- 
spectral density function and the better known cross-correlation function are Fourier 
transform pairs. The primary advantage of using the cross-spectral density function 
is that the correlation can be expressed directly as a function of frequency. The 
mathematics involved in deriving the cross-spectral density function by taking the 
Fourier transform of the cross-correlation function is well documented and will not 
be repeated here. It suffices to state that the resulting cross- spectra I density function 
is a complex quantity composed of real and imaginary parts. The real component, 
denoted as the co-spectral density, is derived by solving the complex function for a 
time displacement of t =0. Thus, the co-spectral density represents a description of 
the in-phase components of the fluctuating pressure field and is equivalent to the 
narrow-band spatial correlation function. The imaginary component, denoted as the 
quad-spectral density, is derived by solving the complex function for a time displace- 
ment of t = l/4f and represents a description of the out-of-phase components of the 
fluctuating pressure field. 

The microphone measurements for locations along the test panel longitudinal center- 
line were analyzed for cross-spectral density characteristics. The method of analysis 
consisted of digitizing successive pairs of microphones in both the upstream and wake 
regions;. For the present analysis, no attempt was made to perform cross-spectral 
density analyses for varying separation distance between microphones except for those 
where the separation between neighboring microphones as installed in the test panel 
varied. Since the majority of the microphones along the longitudinal centerline were 
separated by one inch, most of the analyses were performed for the one-inch separa- 
tion distance. The microphones for the most upstream pair were separated by two 
inches, and the microphones for the two most downstream pairs were separated by 
two inches and four inches. Thus, these microphone pairs are the only ones for which 
cross-spectral densities were computed with separation distances greater than one 
inch. The purpose of this approach was to obtain maximum correlation between the 
fluctuating pressure measurements. Further, the overall fluctuating pressure measure- 
ments indicated that the protuberance flow field was significantly non-homogeneous 
in space, particularly for the upstream separated flow region, and it was felt that the 
correlation lengths would be relatively short. 

Cross-spectral density results for the 8-inch diameter protuberance are presented as 
normalized co- and quad-spectra in Figures 63-68. These results represent areas of 
maximum correlation between neighboring microphones located in the upstream 
separated flow and downstream wake regions for free-stream Mach numbers of 1 .2, 

1 .4 and I .6 and protuberance heights of 1 .0 and 2.0 diameters. It will be noted 
that the region of maximum coherence of the fluctuating pressures for the protu- 
berance induced separated flow region (Figures 63-65) occurred between the point of 
separation and the point of maximum fluctuating pressure, Further, the trend in the 
quad-spectral density and the negative values of the normalized frequency parameter 
show the presence of negative convection velocities, thus indicating that the dis- 
turbances are convected counter to the direction of the free stream. In the wake. 


51 



maximum coherence was observed in the neck and far wake regions. Results presented 
in Figures 66-68 represent the neck region since, here, the separation between micro- 
phones was a minimum of one inch. In the far wake, microphones which were paired 
for the analysis were separated by 2 and 4 inches and some loss of coherence was 
observed at the lower frequencies as discussed later. In general, the wake was 
characterized by positive convection velocities as exhibited by positive values of the 
normalized frequency parameter and the trend in the quad-spectral density. 

A comparison of the various cross-spectra results shows a definite dependency on both 
free-stream Mach number and protuberance height-to-diameter ratio. In general, the 
larger values of and h/D exhibited the greater coherence of the fluctuating 

pressure environment. 

It will be noted that results are not presented for M^ < 1 .0. For this Mach number 

range, the upstream separated flow region was confined to a region in close proximity 
to the protuberance, and the 1-inch separation of the microphones appeared to be 
small relative to the structure of the fluctuating pressure mechanisms. This will become 
clear from the discussion of the model of the separated flow field which is presented in 
Section 6.4. In the wake region, the near wake appeared to extend further down- 
stream for Mqq < 1 .0, than for Mqq > 1 .0. As a result, the region of maximum 
coherence was near or downstream of the most aft microphones and the cross-spectra 
were poorly defined. Thus, the 4-in. diameter protuberance results will be used to 
analyze the wake environments for M^ < 1 .0; however, these data are not available 
for the present report. 

Narrow-Band Convection Velocity 

Before discussing the narrow -band convection velocity results, distinction is made 
between a fluctuating pressure field which is composed of vorticity (or turbulence) 
mode disturbances and that which is composed of acoustic (or sound) mode disturb- 
ances. Vorticity mode disturbances result from a turbulent flow structure which is 
convected at some fraction of the local mean flow velocity. Further, the convection 
velocity is in the same direction as the mean velocity of the turbulent environment. 

On the other hand, acoustic mode disturbances consist of sound waves which piopagate 
in the fluid medium. The propagation velocity relative to the fluid is at the local 
speed of sound with a direction independent of the local velocity. Propagation in the 
direction of the mean flow will have a propagation velocity of the speed of sound plus 
the mean speed; whereas, propagation counter to the flow direction will have a propa- 
gation velocity equal to the difference between the local speed of sound and the local 
velocity. Obviously, acoustic mode disturbances cannot propagate upstream in a 
supersonic flow. 

For the present computation of the narrow-band convection velocities, use was made 
of the narrow band phase angle, as shown in Equation 4.27. Thus, the 

resulting velocities are in actuality the narrow band phase velocities which are the 
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same as the true convection velocities only if the turbulence structure is frozen in the 
moving frame of reference. The idea of a velocity of convection of turbulent fluctua- 
tions stems originally from Taylor's work on grid turbulence (Reference 36). Taylor 
pointed out that, if the turbulence level was low, the time variation in the velocity u 
observed at a fixed point in the flow would be approximately the same as those due to 
the convection of an unchanging, frozen, spatial pattern past a fixed point with the 
mean flow velocity U, i.e., that u(x,t) « u(x - Ut, 0) where x and t represent 
distance measured downstream in the mean flow direction and time, respectively. 

This is known as Taylor's hypothesis. Lin (Reference 37) has shown that this hypothesis 
is valid only if the turbulence level is low, viscous forces are negligible, and the 
mean shear is small . It is unreasonable to expect the hypothesis to apply to the high 
intensity perturbed flow field induced by the protuberance. For the present compu- 
tation, measurements fora constant separation, £, was used and the results represent 
the broad band in wave number. A discussion on convection velocities using wave 
number spectra is given in Reference 38. Due to the non-homogeneous nature of the 
flow field, computations of convection velocity for varying 4 did not seem practical, 
particularly for the upstream separated flow field. 

The variations of the normalized narrow band convection velocity, U (f)/U , with 

c ^ 

normalized frequency, f l/U^, are presented in Figures 69-74. These data corre- 
spond to the locations in the perturbed flow field and the test conditions for which 
cross-spectra results are presented. Also presented in Figures 69-74 are the esti- 
mated maximum local velocities which were computed based on free-stream speed 
of sound, free-stream total pressure and local static pressure. Pitot pressure measure- 
ments were not taken in the perturbed flow field near the protuberance, and conse- 
quently, the values of local convection velocity are only crude estimates of the 
maximum values. These estimates are conservative in that the local total pressure 
was probably significantly less than the free-stream total pressure and thus would give 
a lower local velocity. 

For the upstream separated flow region (Figures 69-71) the narrow band convection 
velocities were negative indicating that the disturbances were generally convected 
(or propagated) counter to the direction of the free-stream. At low values of the 
frequency parameter (low frequencies since I was held constant for the present 
analysis), the convection velocity was a small fraction of the free-stream velocity. 
However, the convection velocity approaches the free-stream value at high fre- 
quencies, and in certain cases, exceeded the free-stream velocity. In every case, 
the estimated local velocity was exceeded at the high frequencies. This suggests 
that appreciable non-frozen turbulence and/or acoustic mode disturbances were 
present in the perturbed flow field upstream of the protuberances at the high fre- 
quencies. The broad -band convection velocity, indicated by the broken line in 
each figure, was either less than or approximately equal to the estimated maximum 
local velocity. It will be noted that the broad band convection velocities are 
heavily weighted to the high frequencies. This characteristic results from the broader 
bandwidths associated with the high frequency data points (see Equation 4.23). It 
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should be noted that the convection velocity data do not reflect the power spectral 
density of the fluctuating pressure energy. Thus, a comparison of the power spectra 
with the narrow band convection velocities will reveal that most of the fluctuating 
pressure energy is convected at a small fraction of the free-stream velocity in a 
direction counter to the free stream. Thus, the true convection velocity from the 
standpoint of convected energy corresponds to the low frequency values. 

In the wake, the narrow-band convection velocities were generally positive indicating 
that the disturbances are convected (or propagated) in the same direction as the free 
stream. Appreciable acoustic mode disturbances appear to be present in the mid- and 
high-frequency ranges as indicated by the large peaks in the narrow -band convection 
velocity curves (Figures 72-74). For certain frequency bands, negative convection 
velocities were computed from the fluctuating pressure measurements which further 
suggests the presence of acoustic mode disturbances in the wake. These data are not 
plotted; however, the frequencies at which they occurred are indicated by broken 
lines between neighboring data points in Figures 72-74. 

6 .4 Analysis of the Protuberance Induced Flow Field 

Most of the technical effort underlying this report has been concentrated in develop- 
ing data reduction techniques and computer programs, reducing the data to tabulation 
and plot form, and evaluating the accuracy of the final values. For this reason, a 
comprehensive analysis of the results in regard to the development of mathematical 
models has not been made. At the present time, additional test results are being 
reduced which will provide more insight into the three-dimensional protuberance 
flow field, and it is felt that these data are required before a comprehensive analysis 
of the flow can be made - particularly in regard to the specification of empirical pre- 
diction formula. For this reason, only a cursory analysis of the protuberance induced 
flow field is presented in this section. 

A composite schematic showing typical axial distributions of static pressure coefficient 
and root-mean-square fluctuating pressure coefficient, and a possible model of the flow 
field upstream of the protuberance is presented in Figure 75. Sufficiently upstream of 
the protuberance, the flow is attached which results in static pressure coefficients of 
approximately zero, and fluctuating pressure coefficients of approximately 0.01 . The 
latter results primarily from tunnel background noise as discussed in Section 5.0. For 
the case shown, which is the = 1 .60, h/ D = 2 .0 condition, the most upstream 

boundary of the perturbed flow is the shock wave. The shock wave induces a nonlinear 
increase in static pressure coefficient and a relatively small peak in the fluctuating 
pressure coefficient. The mean shock wave location corresponds to the inflection 
point in the static pressure data. Immediately aft of the shock wave, the boundary 
layer undergoes separation which, in turn, results in a region of relatively constant, 
plateau, static pressure. In the plateau region, a gradual increase in fluctuating 
pressure level with decreasing distance from the protuberance is observed. In con- 
trast to these data, the two-dimensional protuberance (or step) case shows a region 
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of plateau fluctuating pressure corresponding to the region of plateau static pressure 
(see References 27 and 35 for examples) . This latter discrepancy between two- and 
three-dimensional separated flows is the first indication (as the flow is examined 
moving downstream toward the protuberance) of dissimilarities between the two cases. 
Further, the fluctuating pressure levels for the three-dimensional case are an order 
of magnitude greater than that for the two-dimensional separated flow region. Im- 
mediately downstream of the plateau region, the flow exhibits a sudden reduction in 
static pressure, and the beginning of this sudden reduction corresponds to a peak in 
the fluctuating pressure level . Corresponding to the region of low static pressure is, 
what may be envisioned as a plateau region of fluctuating pressure for the M^ = 1 .60 

condition; however, this is not typical of all Mach numbers. For lower Mach numbers, 
the fluctuating pressure decreases toward the protuberance for the region of low static 
pressure as discussed in Section 6.3.2. The distinct nature of the sudden static pres- 
sure reduction and the strong levels of fluctuating pressure are unique to the three- 
dimensional separated flow case. Before describing the model of the separated flow 
region, some additional features of the flow field as given by oil flow patterns of the 
separated region will be discussed. 

A typical oil flow picture. Figure 26, was discussed briefly in Section 6.3. Colored 
oil was injected into the attached boundary layer upstream of the perturbed flow 
region as well as in close proximity to the protuberance. The oil from the three 
upstream orifices flowed aft and was deflected around the protuberance by the sepa- 
rated flow field. The oil injected into the separated flow region, flowed in a direc- 
tion counter to the free stream - indicating a reverse flow region, as well as transverse 
to the free stream - indicating side relief due to the three-dimensional nature of the 
flow field. The important feature of the oil flow pattern in the separated flow region 
is the concentration of oil along a secondary front which is downstream of the front 
caused by separation of the boundary layer. This suggests a condition of low velocity 
at angles perpendicular to the secondary front. The presence of this low velocity 
region imbedded within the separated flow has not been observed in two-dimensional 
separated flow and is apparently unique to the three-dimensional case. The oil flow 
pictures also reveal the flow patterns in the wake region; however, an adequate dis- 
cussion of this region was presented in Section 6.3. 

Based on the characteristics of the induced separated flow as exhibited by the static 
pressures, fluctuating pressures, and oil flow pictures, a model of the upstream 
perturbed flow can be postulated. Consider a system of horse-shoe vortices, which, 
when viewed in the plane of symmetry appear as shown in Figure 75. Two major 
vortices are assumed to be rotating such that a condition of reverse flow exists near 
the wall of the test panel with a third, smaller, vortex trapped near the wall between 
the two major vortices and rotating in the opposite direction. Thus, the boundary 
conditions at the inner region of the separated shear layer and within the separated 
flow are satisfied. Also observe that the upstream vortex may be elongated which 
would result in a relatively constant static pressure, plateau, region. The aft vortex, 
because of the geometry of the flow, may be more circular which would induce a 
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nonlinear variation of the wall pressures as observed in the static pressure distribu- 
tion. Further, because of the different static pressure levels associated with the two 
vortices, it may be postulated that the inner vortex has a higher tangential velocity 
at the wall than does the forward, elongated vortex. Finally, the shearing mechanism 
that would result due to the opposite velocity vectors in the region between the two 
vortices would cause a discontinuity in the wall shear stress which would enable oil 
to be trapped in this region (corresponding to the secondary front noted in the oil flow 
pictures) and also noise resulting from the shear interaction could generate high level 
fluctuating pressures as noted in the axial distribution of ACp(RMS), Figure 75. The 
direction of rotation of the vortices is verified from the oil flow patterns over the 
surface of the test panel (Figure 26) as well as in photographs of titanium dioxide 
studies of the protuberance flow pattern. Figure 76. It is interesting to note that the 
titanium dioxide flow patterns (side view) shows the separation of the flow from the 
wall of the protuberance near the 90-degree ray very clearly as was observed in the 
static and fluctuating pressure measurements over the wall of the 8-inch diameter 
protuberance . 

Other data which substantiates the assumption of a double vortex system are the 1/3- 
octave spectrums. If was pointed out that, based on the static pressure levels, the 
aft vortex appears to have the greater rotational velocity. Thus, fluctuating pressures 
being convected past the microphones in the region of the aft vortex would be ex- 
pected to have energy concentrated at higher frequencies than those associated with 
the upstream vortex. Typical 1/3-octave spectra show this trend as shown in Figure 77. 
Further, cross-spectral densities for microphone pairs located beneath the separated 
flow were obtained at various locations in the separated flow regions and these data 
yield even greater insight into the mechanism of the flow. If the postulated model is 
accurate, good coherence should be shown for microphone pairs located in either of 
the two regions, i .e., microphone pairs located either upstream or downstrem of the 
peak ACp(RMS). However, poor coherence should be evident from cross-spectral 
densities for microphone pairs which are separated by the peak since the microphones 
would sense fluctuating pressures from two separate vortex systems. 

Cross-spectra for a representative condition (Mqq = 1 .60, h/ D = 2.0) are presented 

in Figure 78. It will be noted that reasonably good correlation is shown for micro- 
phone pairs located beneath the forward vortex, poor coherence for microphone pairs 
located near the region between the two major vortices, and fair correlation for micro- 
phone pairs located in the region of the inner vortex. It will be noted that distur- 
bances generated in the region of interaction between the two major vortices would 
not be convected directly across the microphones near the protuberance because of 
the direction of rotation of the inner vortex. Thus, some loss of coherence (compared 
to that for the outer vortex system) may be expected. 

One final sample of data which adds support to the validity of the postulated model 
of the upstream separated flow is presented in Figure 79. Here, the flow pattern 
upstream of a cylindrical protuberance has been photographed using thin filaments 
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of smoke for flow visualization (from Reference 39). The multiple vortex pattern 
upstream of the protuberance is clearly defined. This condition was extremely low 
speed (on the order of 5 ft/sec); however, the present data suggests that the basic 
features of the flow field remain unchanged at transonic Mach numbers. 

The flow structure in the wake consists of a converging flow region, a neck region of 
nearly parallel flow, and a diverging flow region when viewed in the plane of the 
test panel . This flow field is further complicated by the flow over the top of the pro- 
tuberance, particularly for h/D <1.0. The characteristics of the fluctuating pressure 
field have been previously discussed; however, some additional insight can be gained 
from the cross-spectra for various locations in the wake. Figure 80. In the near wake, 
r/ D < 1 .00, the cross-spectra decay verv rapidly and in general show poor correlation. 
The correlation increases to a maximum in the neck region and this is attributed to the 
non-dispersive nature of the parallel flow in this region. In the far wake, the flow has 
a diverging pattern; however, good correlation was still attained. The larger separa- 
tion of the microphones in the far wake resulted in a loss of coherence in the low 
frequency components. 

In summary, the preliminary analysis of the test results have revealed some important 
characteristics of the static- and fluctuating-pressure environments induced by three- 
dimensional protuberances. Based on this analysis, some conclusions have been drawn 
in regard to the perturbed flow field and areas for further, expanded, studies have been 
defined. These points are discussed in the next section. 


7.0 


CONCLUDING REMARKS 


7.1 Conclusions 

A wind tunnel investigation has been conducted to study the flow field induced by 
three-dimensional protuberances at transonic Mach numbers. Selected results for 
the generalized protuberance - which consisted of right circular cylinders - have 
been presented and discussed. In particular, the static pressure field, the fluctuating 
pressure field, and general features of the perturbed flow have been evaluated. From 
this study, the following conclusions have been made. 

7.1 .1 General Characteristics of the Protuberance Induced Flow Field 

1 . The protuberance flow field consisted of two predominant regions - an up- 
stream separated flow region and a downstream wake region. The upstream 
separated flow region was composed of a complex, multiple, vortex flow 
structure. This flow structure appears to be unique to the three-dimensional 
separated flow case and has little similarity to the two-dimensional, step 
induced, flow structure. A postulated model of the three-dimensional 
separated flow region consists of two major vortices rotating in the same 
direction with velocity vectors at the wall which are counter to the free- 
stream direction. Static- and fluctuating-pressure measurements and oil- 
flow patterns appear to agree with the postulated flow structure. The wake 
region consisted of a near wake of converging flow, a neck region of 
approximately parallel flow , and a far wake of diverging flow. For small 
ratios of protuberance height to diameter, the wake flow structure was com- 
plicated by the flow over the top of the protuberance which appeared to 
re-attach on the surface supporting the protuberance in the vicinity of the 
neck of the wake . 

2. Increasing either protuberance height or free-stream Mach number caused an 
increase in the upstream separation length and a corresponding growth of the 
separated flow structure. 

3. The perturbed flow field scales approximately with protuberance diameter, 
although some trends in the normalized data with protuberance size and 
Reynolds number were detected in the static- and fluctuating-pressure 
measurements . 

7.1 .2 Protuberance Induced Static Pressure Field 

1 . The protuberance induced static pressure field upstream of the protuberances 
exhibited the following characteristics: 

a. There was an increase in static pressure at the outer extremities of 
the upstream perturbed flow field that induced boundary layer 
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separation. At M Q0 < 1 .20, the pressure gradient was relatively 
small; whereas, at M^^ 1 .20, the pressure gradient was large and 
resulted from the pressure discontinuity of the shock wave. 

b. The initial pressure gradient was followed by a region of nearly 
constant, plateau, static pressure which corresponded to the outer 
region of the separated flow field. 

c. A region of low static pressure (below the plateau level) was 
exhibited by the inner region of the separated flow field. 

d. At the base of the protuberance, maximum static pressure was 
reached . 

e. Variations in protuberance height and free-stream Mach number 
caused similar variations in the region exposed to the protuberance 
induced static pressure field. An increase in the test variables 
caused an increase in the area of the protuberance induced static 
pressure field. 

2. The protuberance induced negative static pressure coefficients immediately 
downstream of the protuberance which corresponded to the near wake and 
neck regions of the wake flow field. Approximately free-stream static 
pressure was measured in the far wake. 

3. The static pressure field over the wall of the protuberance exhibited 
noticeable variations in both the vertical and circumferential directions 
over the forward half of the cylinder. Variations in the vertical direction 
corresponded to the various flow regimes of the inflowing perturbed flow 
field. Circumferential variations indicated that the flow separated from 
the wall of the protuberance near the 90-degree ray. 

4. Variations in Reynolds number and protuberance size caused small variations 
in the wake pressures at subsonic Mach numbers, and small, but consistent, 
variations in the upstream pressures at supersonic Mach numbers. 

7.1 .3 Protuberance Induced Fluctuating Pressure Field 

1 . The protuberance induced fluctuating pressure field upstream of the protu- 
berances exhibited the following characteristics: 

a. There was a gradual increase in fluctuating pressure in the axial 

direction which started at the outer extremity of the perturbed flow 
field and reached a peak near the mid-point of the separated flow 
region. The peak levels of fluctuating pressure corresponded to a 
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region between two major vortices and appeared to result from the 
shear interaction of the vortices. Peak levels were in the range 
from 12 to 17 percent of the free stream dynamic pressure and are 
an order of magnitude greater than that observed for two-dimensional 
separated flows. 

b. The inner region of the separated flow field exhibited a high level of 
fluctuating pressure which corresponded to the inner vortex system. 
For M^ = 1 .60, the fluctuating pressure levels for the inner region 

were approximately constant and equal to the peak level in the 
separated flow field. For M^< 1 .60 , the fluctuating pressure levels 

in the inner region decreased below the peak level . 

c. For = 1 .40 and 1 .60, shock -boundary interaction upstream of 

the separation point resulted in a small peak in the fluctuating pres- 
sure levels for the region beneath the shock wave. 

d. Variations in protuberance height and free-stream Mach number 
caused similar variations in the extent of the separated flow field, 

i .e., increasing test variables caused an increase in the size of the 
separated flow field and, thus, the region exposed to protuberance 
induced fluctuating pressures. Otherwise, the separated flow struc- 
ture and the corresponding fluctuating pressure field appeared to 
remain relatively invariant with changes in protuberance height and 
free-stream Mach number. 

e. The power spectra corresponding to the peak levels of fluctuating 
pressure exhibited a high energy level at low frequencies with a 
roll-off at high frequencies of approximately 6 dB per octave. 

f. Maximum correlation of the fluctuating pressures within the separated 
flow field corresponded to the region of the outer vortex system. 

g. Fluctuating pressures within the separated flow field were convected 
in a direction counter to that of the free stream. 

h. The fluctuating pressure field was largely non-homogeneous and this 
is attributed to the complex vortex structure of the separated flow 
field . 

2. Fluctuating pressures over the wall of the protuberance exhibited peak values 
in the region of separation of the flow from the wall of the protuberance and 
in the wake region of the protuberance at subsonic Mach numbers. At super- 
sonic Mach numbers, maximum levels of fluctuating pressure occurred over 
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the front half of the protuberance corresponding to the region of impingement 
of the upstream separated flow field. At = 0.60, the flow separated from 

the wall of the protuberance at the 90-degree ray and resulting fluctuating 
pressures exceeded 170 dB. 

3. The protuberance induced fluctuating pressure field downstream of the protu- 
berances exhibited the following characteristics: 

There was a moderate increase in fluctuating pressure levels in the 
axial direction within the near wake, peak levels of fluctuating 
pressure in the neck region and gradually decreasing fluctuating 
pressure levels in the far wake. 

Power spectra of the fluctuating pressures revealed discrete fre- 
quency energy at certain protuberance heights and Mach numbers. 

The Strouhal number, based on protuberance diameter, associated 
with the peaks in the power spectra was approximately 0.28. 

Maximum correlation of the fluctuating pressures within the wake 
flow field corresponded to the neck and far wake regions. 

The fluctuating pressure characteristics of the wake appeared to vary 
with both protuberance height and Mach number . 

7.2 Areas for Further Study 

Based on the foregoing analysis of the test results, important areas requiring further 

analysis can be defined. These areas are as follows: 

1 . A more extensive study of the spectra and cross spectra of the fluctuating 
pressures recorded by the centerline distribution of microphones. These 
data will be analyzed to define the structure of both the upstream separated 
flow (i .e., verification or disproval of the postulated model of the separated 
flow) and the wake flow aft of the protuberance. 

2. Examination of the microphone and static pressure data for locations around 
the protuberance in the plane of the test panel . These data will be analyzed 
to define the generating mechanisms associated with the three-dimensional 
flow around the protuberances that caused differences between two- and 
three-dimensional flow patterns as observed in the centerline axis data. It 
is felt that the transverse velocity and acceleration components induced by 
the three-dimensional protuberance shape are the key to the multiple vortex 
flow pattern which has been postulated. These data will also provide a 
basis for the development of prediction techniques for the steady and fluctua- 
ting air loads around the protuberances. 
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Further examination of the microphone and static pressure data recorded for 
sensors distributed over the wail of the 8-inch diameter protuberances. These f 

data will be analyzed to define the steady and fluctuating air loads on pro- f 

tuberances and to assist in defining the three-dimensional flow patterns. 

4. A more detailed study of the effects of Reynolds number and protuberance 

scale on the flow field particularly for the fluctuating pressure environment. f 

5. An evaluation of the free-interaction hypothesis as it applies to three- 
dimensional flows. This hypothesis, advanced by Chapman, Kuehn and 
Larson, Reference 2, suggests that turbulent boundary layers undergo separa- 
tion in a manner independent of the detailed geometry of the cause. Al- 
though proven for two-dimensional separated flows, its application to three- 
dimensional flows has not been verified. If the hypothesis is true for the 
three-dimensional case, prediction of protuberance flow fields for specific 
geometries can be made based on data obtained for generalized protuberance 
shapes. Further, this hypothesis applies only to the upstream separated flow 
region. It is also of interest to determine if the free-interaction concept 

can be applied to the wake region as well . 'X 

In addition to those described, it is anticipated that other areas requiring study will 
become evident during the course of the analysis. 
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APPENDIX A 


DATA REDUCTION EQUATIONS FOR MEAN FLOW PARAMETERS 


The parameters used in the specification of the mean flow characteristics are presented in terms 
of the measured static and total pressures as follows. 


1 * Local Static Pressure Coefficient 


2 . 


c p= q 


“ -'co 


(A— 1) 


oo 


where Pg and P^ are the local and free-stream static pressures, respectively, and 


q^ is the free stream dynamic pressure. 


Boundary Layer Profile Mach Number (Computed at Both Forward and Aft Rake 
Posi ti ons) 


For P g /Pf n < 0.528, the Rayleigh Pitot Equation was solved for Mach numbers as 


follows: 
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where P $ is the test panel static pressure nearest the rake and Pj is the pitot pressure 
measured by the rake. 

For P s /Pt > 0.528, the standard compressible flow equation was solved, i.e.. 
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3. Boundary Layer Profile Velocity Ratio (Computed at Both Forward and Aft Rake Positions) 
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where U s and M $ are the velocity and Mach number corresponding to the local static 
pressure P , and U n and M n are the profile velocities and Mach numbers corresponding 
to the measured pitot pressures, Pj. . 

4. Boundary Layer Displacement Thickness (Computed at Both Foreword and Aft Rake 
Positions) 



(A-5) 


where 5 is the y position for the pitot probe that is the furtherest from the test panel for 
each rake, p^ and p^ are densities corresponding to P $ and P^. 


5. Boundary Layer Momentum Thickness (Computed at Both Forward and Aft Rake 
Positions) 
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6. Empirical Boundary Layer, Displacement and Momentum Thicknesses (Computed at 
Model Station Zero) 


where 


5 * = 


X [1.3 + 0.43 M^] [5/X] 


10.4 + 0.5 M^ [1 + (2) (10" 8 ) Re^-X] 1 ^ 


oo 


6/X = 0.37 f Re • X1+ 5 

OO 


1 + 


Re -X 
00 

6.9 • 10 7 


Vio 


and 


X = 124.36 inches 


0 = 


1.3 + 0.43 M^ 


(A-7) 


(A-8) 


(A-9) 


68 


test fart point t i he cate *edc PROPULSION wind tunnel 

Tf-184 * 154 5 251 6- 7-68 ' " ~ ' TRANSONIC 1*>T 

Ml PTT PI Cl 2F RXto-6 NC* P “T A 'PTb fcJ PC8 PTT wo TPR w r/1wT^f^ fc “'T^T'TT V”TT’J 

1.5985 1572-1 37r>.7 663.0 41227 3.000 45 1572-1 1572.3 3P4.5 384.3 1060.0 225.2 1.4P3 0.0480 0.51 0.51 117 104 

K1 K2 K3 M 55 5J K7 *8 K9 *10 CONF 


0.2083 0 

O 

♦ 

© 

• 

.2120 

C • 2094 

0.2094 

0.203 

5 0.0000 0.0000 0 

.0000 0-2037 

06 





PHI 

0 

15 

30 

Tb0-INCP PRCTUbtRANC6 X0D= 2.00 CORRECTED 

45 67.5 90 112.5 135 150 157.5 165 

170 

175 

180 

270 

R 

1.75 

1.531 

0.332 

1.583 

0.335 

1.589 
C -236 

1 .596 
0.341 

1.512 

0-296 

1.306 
0.175 - 

0-906 

0-053 

0-591 

-0-234 

0.507 0.483 0-437 

-0-237 -0.300 -0.329 - 

0.421 
0.334 - 

0.572 C.630 

0.250 -O .216 

- 

2.75 

1.253 

0.146 

1.229 

0.130 

1-177 

0-101 

1.176 

0.102 

1.2n9 

9.119 

1.232 

0.133 

1-040 

0.023 

0.697 

-0-176 

3.547 0.494 

-0.262 -0.293 

0.449 

-C.320 - 

0.500 
0.282 - 

0.647 

0.202 

0.73^ 

-0.153 


3.75 

1.504 

0-282 

1 .489 
0.27e 

1 .457 
0.261 

1.397 

0-227 

1.254 

0.145 

1.088 

0.050 

1-061 
C • 034 

0-738 

-0-151 

0.551 0.479 

-0.260 -0.295 

0.486 
-0.294 - 

0.613 
0*219 - 

0-670 

0.188 

0.736 

-0.149 

•• ■ 

v.vr 

1-314 

0.165 

1.466 

0.245 

1.453 
C .239 

1.446 

0.235 

1.355 

0.187 

1 .207 
0.110 

1-036 

0-019 

o.9oT“ 

-0-057 

0.533 0.4b3 

-0-249 -0.290 

0.532 
-0.249 - 

0.6 33 
■0-197 - 

0.697 

•0-160 

0.757 

-0.13? 

1.1*7 

0.101 

5.50 _ 

1-073 
' 0-038 

1.374 

0.194 

1.383 

0*200 

1.389 
0 • 2 G 5 

1.340 

0.180 

1.232 

0.323 

1-031 

0-017 

0-79n 

-0-117 

0.533 0.454 

-0.249 -0.239 

0.674 
-0.199 - 

0-703 
■0.158 ■ 

0-795 

■0-108 

0.840 

-0.090 

1.232 

0.123 

6.50 

1.003 

0*002 

1.137 

0-071 

1.170 
C • 090 

1.270 

O.l-il 

1.370 

0.168 

1.243 

0.128 

1-055 

0-029 

0.806 

-0-10? 

0.549 n.462 

-0.239 -0.287 

0.776 
-0.118 - 

0.974 

■0.066 

1-009 

0-005 

0.996 

-0.002 

1.235 

0.125 

7.50 

1.007 

0.003 

1-01? 

0-006 

1.010 

0.005 

i.oe9 

0.046 

1.241 

0.176 

1.231 

0.127 

1-087 

0-046 

0-804 

-0-10? 

0.540 0-629 

-0.247 -0.197 

0.959 

-0.021 

l.«l“ 

0-009 

1.068 

0-036 

1 .105 
0.059 

1.244 

0.127 

STfcT 

1.004 

0-002 

1.008 

0-004 

1.001 

C-000 

1.014 

0-OC7 

1 • 1 «5 
0-075 

1 .24.1 
r. .125 

i . 1 3<r~ 

0-069 

0 - 3o8 
-0-101 

0.593 0.558 

-0-216 -0-076 

1.018 

0.009 

1 . oTT - 
0.039 

1.0 6^2 
0.0 33 

1.031 1.232 

0.046 0.120 

9.50 

1-006 

0-003 

1-006 

0-002 

C .991 
-G .005 

0 .9G8 
-o-oci 

1.061 

0-027 

1.717 

n.113 

1-145 

0-075 

0-827 

-0 .074 

0.591 0.914 

-0.219 -0-046 

1.028 

0.015 

1 .055 
0-033 

1.070 

0.038 

1.067 

0.033 

1.201 

0.132 

IOTTo 

0.998 1.005 

-0-001 0-002 

1*004 

0-002 

1-010 
0 .005 

1.019 

0*010 

l.j9l 

0.097 

1-148 

0-077 

•.356 

-0-075 

0.67i 0T941 

-0.177 -0-032 

1.044 

0.023 

l.C 42 
o • n 33 

1 • 05T 
0.078 

1.65*“ 

0.032 

i.i7a 

0.066 

11.50 

0.995 

-0-002 

1.006 

0.003 

0 .998 
-C -001 

1 . C 0 3 
0.002 

1.008 

0.004 

1.357 

o.ofiu 

1-163 

0-080 

0.856 

-3.077 

0.61 7 0.953 

-0.093 *0.025 

1.023 

0.012 

1.060 

0.032 

1-061 

0.033 

1.054 
0.0 31 

1.121 

0.361 

f?750 

l.aii 

0-006 

1.006 

0.003 

1-014 
C -007 

0.997 

-0.002 

1.003 
C .002 

1 . P85 

0.044 

1 167 
0-081 

~'.v:6 7i 

-0 . 069 

0.86 3 P.976 

-0.074 -0.013 

1.028 

C.015 

1 . P 37 
0.019 

1.055 

0.029 

1.051 

0.029 

1 .072 
0.036 

13.50 

1.002 

0-001 

1.007 

0.0»*.3 

0.994 

-0-003 

1 .002 
0.001 

0.993 

-0.004 

1.060 

0.031 

1-146 

0-076 

0.876 

-0.065 

0.890 0.981 

-0.059 -O.flilO 

1.01 8 
0.009 

1.039 

0.020 

1.047 

0.025 

1.041 

0 .023 

1.078 

0.039 

14 . 5 b 

1.000 

- 0-000 

1.044 

0.022 

0.993 

-0-004 

0.997 

-0.00? 

1.006 

0.003 

1.037 

0.019 

1-149 

0-076 

0.39* 

-0.057 

0.9l4 0.994 

-0.046 -0-003 

1.024 

0.013 

1.024 

0.013 

1.036 

0.019 

1 .04n 
0.02? 

1.074 

0.012 

15.50 

1.003 

0-001 

1.144 

0.072 

t .985 
-0 -007 

1.006 
0 • 0 0, 3 

1.008 

0-004 

0.992 

-0.004 

1-139 
C -069 

0.904 

-0-051 

0.926 1.003 

-0.039 P.002 

1.026 

0.013 

1 .026 
0.013 

1.040 

0.021 

1 .04? 
0 . 02 T 

1.077 

0.013 

16.50 

1.006 

0-003 

1 .477 
0.249 

C .978 
-€ -Oil 

0.995 

-0.002 

0.993 

-0.004 

0.996 

-0.092 

1-130 

0-065 

0-9n« 

-0.049 

0.939“ n.976 

-0.032 -0-013 

1.0 ?•< 
0.012 

1.019 

0.010 

1.013 

0.007 

1.024 

0.013 

1.001 

o.ono 


TEST 

184 


part 

154 


POINT 

5 



R = 

1775b 

18.50 

19.50 

2 0.50 

21 . 5o 

~2?75o~ 

23759 — 

24.50 25. 50 25750 ?3.00 




PHI «0 

1.008 

1.005 

C .987 

C .999 

1.004 

1.006 

0.999 

6.995 0.999 0.995 ' 






; 

0.003 

0.002 

-0.006 

-C.OOC 

0.002 

C.003 - 

0.001 - 

0.003 -0.001 -0.002 




* 

PHI *180 

1.020 

1.016 

1 .015 

1 .021 

1.016 

1.024 

1.019 

1.019- 1.024 1.023 1.025 






0.011 

0-009 

0-008 

C . C12 

0.009 

0-013 

0.011 

6. oil 0 . 01 3 o.oiJ 6.014 




' PIN 

HP 

" 12P 13P 

14P 

15P 

••••COMPARISON 

1 PRESSURES**** " iftp "■ 17P ' ~ : 

l*P 

19P 

10P 

530.5 

544. 

5 541.6 5 

43.9 

54J.9 

542.5 



528.0 0.0 

0.0 

0.0 

528.3 

PTOT 

PS1 

PS2 










368.2 

371. 

4 374.9 










’ -- - 

- . _ 

- - - * * 

- 


— 

"FIXED 

RAKE ( 2 > 

* ~ 

. 

...... -- 

------ - 

Ml 

M2 

M3 

M4 

N5 

M6 

H? 

M9 

M9 

Mto Mil M 12 M13 Ml 4 

M15 



1.000 

1.213 

1.265 

1.299 

1.317 1.333 

1.347 

1.372 

1-400 

t-42? 1.4*7 1.464 1.478 1.494 

1.508 



H16 

H17 

818 

MIS 

"20 

"21 

*22 

M23 

824 

M25 




„ _1-519 

1.529 

1.539 

1.539 

1.53; 

e 1.531 

1.537 

1 .543 

1.548 

1.55? 

, - . — - - - 

... - 

. - 

Ul/U 

U2/U 

U3/U 

U48L 

U5/U 

U6/L 

U7/U 

U8/U 

U9/IJ 

U10/U Ull/U U12/U U1 3/U U14/U 

U15/U 



0.702 

0.820 

0.8*7 

0-86 4 

0.87; 

2 c.aei 

0.887 

0.899 

0-912 

0.923 0.934 0.942 0.948 0.95* 

0.961 




U16/U U17/U 018/U UlS/U U20/1 U21/U U22/U U23/U U24/U IJ25/H 

0.966 p.«7l 0-975 0.975 0.97? O^"?? 0.974 0-977 0-979 0.981 


DISPLACEMENT THICKNESSES 
12 3 

0 . 0 000 0 . 2892 0.2 957 


MOMENTUM THICKNESSES 
1 2 3 

0.0000 0.1981 n.1245 





jfi 



APPENDIX B 


DATA REDUCTION EQUATIONS FOR UNSTEADY FLOW PARAMETERS 


The statistical functions used in the specification of the fluctuating pressure characteristics of 
the unsteady flow field are presented in analog form. It is felt that the analog form of the 
descriptive statistical parameters will facilitate an understanding of the physical significance 
of the functions more so than the digital form. A comprehensive presentation of the digital 
equations used in the present digital computations are given in References 23 and 24. A brief 
discussion of the equations is presented in Section 4.4. The analog forms are as follows. 


1 . Overall RMS Fluctuating Pressure 
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P 2 (t) dt 
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where P(t) denotes the time history record of pressure. 
2. Power Spectral Density 
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where P(t / f / Af) is the filtered portion of P(t) in the frequency range of f - Af/2 to 
f + Af/2 . 


3. Cross Spectral Density 


G(£,f) = C(M) - « Q(£,0 (B— 3) 

where C(£,f) and Q(|,f) denote the co-spectral density and quad-spectral density, 
respectively. The co- and quad-spectral densities are defined as follows: 
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where P(t,f, Af) and P^ (t,f, Af) pressure time histories taken at points separated in 
distance by £ and P|* (t,f,Af) denotes a 90-degree phase shift from P^(t,f,Af). The 
normalized co- and quad-spectral densities are given by: 

C(l,f) 
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and 


4. Phase Angle 


5. Coherence 
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6. Narrow-Band Convection Velocity 


u c (f) = 
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Figure 1 . Schematic of the Test Article Installation in the AEDC - 16T Wind Tunnel Facility 
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Figure 2. Photograph of the Model Installation in the AEDC-16T Wind Tunnel Facility 









Figure 4. Photograph of the Protuberance Drive System and Instrumentation 
Supported Beneath the Test Panel 
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Figure 5. Details of the Cylindrical Protuberances 
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Side View Section 


- Figure 6. Details of the Cylindrical Protuberance Drive System 
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re 7 . Concluded 


(b) Insert Plate for 2- Inch Diameter Protuberance 


(a) Insert Plate for 4- Inch Diameter Protuberance 


Figure 8 . Instrumentation 


Locations in the Insert Panels 
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with Extension 


Figure 9. Instrumentation Location in the 8-Inch 
Diameter Protuberances 
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Figure 10 * Defails of the Boundary Layer Profile Rakes 
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Figure 12, Block Diagram of Acoustic Recording System 













b. Dynamic Pressure 


Figure 13. Variations of Unit Reynolds Number and Dynamic Pressure With 
Free-Stream Mach Number 
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Figure 15. Characteristics of the RMS Fluctuating Pressure Coefficient for the 
Undisturbed Flow Environments Showing the Contributing Sources 












Frequency, f, Hz 


a . M =0 .60 
co 

Figure 16. Characteristics of the Power Spectra for the Undisturbed Flow Environments 
Showing the Contributing Sources 
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Figure 16. Continued 
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Frequency Response Function, f:(f) 



Figure 17. Frequency Response Characteristics of the Microphone/ Data Acquisition System Electronics 














Figure 20= Comparison of the Analog and Digital Computations of 
RMS Fluctuating Pressure Level 
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Figure 21 . 


Variation of Statistical Accuracy Parameters with Frequency for the 
Present Digital Computation of Power Spectral Densities 
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Figure 22 » 


Static Pressure Distributions Along the Longitudinal and Transverse 
Centerlines for Undisturbed Flow over the Test Panel, Re/it = 3x 10 6 
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Figure 24. Variation of Boundary Displacement Thickness 
With Mach Number for Undisturbed Flow 
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Figure 27. Variation of Separation Length with Protuberance Height 
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Figure 28, Variation of Static Pressure Coefficient Along the Longitudinal 
Centerline of the Test Panel, 2-inch Diameter Protuberance, 

Re/ft = 3 x 10 6 
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Figure 28, Continued 
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Figure 28 » Continued 
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Figure 28. Continued 
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Figure 28. Continued 
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Static Pressure Coefficient 



Figure 28. Concluded 
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Static Pressure Coefficient 



Figure 29„ Variation of Static Pressure Coefficient Along the Longitudinal 
Centerline of the Test Panels 4-inch Diameter Protuberance, 
Re/ ft = 3 x 1 0 6 
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Figure 29. Continued 
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Figure 29 . Continued 

20 


1 



-8-6-4 -2 0 2 4 6 8 


0 + = 0° — - 9 + = 180° 

Axial Location, r/D 

e. M = 1.10 
oo 

Figure 29 . Continued 
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Static Pressure Coefficient 



Figure 29. Continued 
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Figure 29 . Continued 
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Figure 30 „ 


Variation of Static Pressure Coefficient Along the Longitudinal 
Centerline of the Test Panel , 8-inch Diameter Protuberance, . 


Re/fi = 3 x 10 6 
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Figure 30. Continued 
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Figure 30= Continued 
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Figure 30. Concluded 
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Figure 31 . Variation of the Separation-to-Free-Stream Pressure Ratio 
with Mach Number 
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Figure 32. Effect of Reynolds Number on the Axial Distributions of Static Pressure 
Coefficient, 8-inch Diameter Protuberance, h/D = 2.0 
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Figure 32. Concluded 
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Figure 33, Effect of Protuberance Size on the Axial Distribution of Static Pressure 
Coefficient, Re^ = 0 „5 x 10 6 , h/D = 2 .0 
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Figure 33 o Concluded 
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0 8-Inch Diameter Protuberance 
O 4-Inch Diameter Protuberance 



Axial Location, r/D 


Figure 34 „ Effect of Protuberance Size on the Axial Distributions of Static Pressure 
Coefficient, Re^ = 1 ,0 x 10 6 , h/D = 2,0 
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Figure 34 „ Concluded 
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c. h/D = 1.5 


d. h/D = 2.0 


Figure 36. Variation of Static Pressure Coefficient Over the Wall of the 8-inch Diameter 

Protuberance,, M = 0.80, Re/ft = 3 x 10 6 
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Figure 38= Variation of Static Pressure Coefficient Over the Wall of the 8-inch Diameter 
Protuberance,. M = 1 = 20, Re/ft = 3 x 10 6 
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Figure 39 „ Variation of Static Pressure Coefficient Over the Wall of the 8-inch Diameter 
Protuberance, M = 1 °40, Re/ft = 3 x 10 6 
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Figure 40. 


Variation of Static Pressure Coefficient Over the Wall of the 8-inch Diameter 
Protuberance, M = 1 .60, Re/ft = 3 x TO 6 
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Figure 41 . Variation of the RMS Fluctuating Pressure Coefficient Along the Longitudinal 
Centerline of the Test Panel, 2-inch Diameter Protuberance, Re/ft = 3xl0 6 
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Figure 41 „ Continued 
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Roof-Mean-Square Fluctuating Pressure Coefficient, AC (RMS) 
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Figure 41 „ Continued 
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Figure 41 . Continued 
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Figure 42o Variation of the RMS Fluctuating Pressure Coefficient Along the Longitudinal 
Centerline of the Test Panel, 4-inch Diameter Protuberance, Re/ft = 3xl0 6 
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Figure 42. Continued 
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Figure 42. Continued 
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Figure 42. Concluded 
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Figure 43. Variation of the RMS Fluctuating Pressure Coefficient Along the Longitudinal 
Centerline of the Test Panel, 8-inch Diameter Protuberance, Re/ft = 3x 10^ 
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Figure 43 o Continued 
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Root- Mean- Square Fluctuating Pressure Coefficient, iC D (RMS) 
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Figure 43. Concluded 
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Figure 44 „ Comparison of the Axial Distributions of RMS Fluctua 
For the Various Cylindrical Protuberances, Re/ft = 3 
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Figure 44 „ Continued 
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Figure 45. Variation of RMS Fluctuating Pressure Coefficient Over the Wall 
of the 8-inch Diameter Protuberance, M = 0.60, Re/ft=3x 10 6 
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Figure 46. Variation of RMS Fluctuating Pressure Coefficients Over the Wall 
of the 8-inch Diameter Protuberance/ M =0.80, Re/ft = 3 x 10^ 
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Figure 47. Variation of RMS Fluctuating Pressure Coefficients Over the Wall 
of the 8-inch Diameter Protuberance, M =1 .00, Re/ft = 3 x 10 6 
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Figure 48. Variation of RMS Fluctuating Pressure Coefficient Over the Wall 

of the 8-inch Diameter Protuberance, M =1 .20, Re/ft = 3 x 10* 
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Figure 49. Variation of RMS Fluctuating Pressure Coefficient Over the Wall 
of the 8-inch Diameter Protuberance, M = 1 .40, Re/ft = 3 x 10 6 
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Figure 52. Power Spectra of the Protuberance Induced Fluctuating Pressure Field 
Upstream (0“ = 0°), 8-inch Diameter Protuberance, M = 0.80, 
Re/ft = 3 x 10 6 00 
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Figure 52. Concluded 
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Figure 53. Power Spectra of the Protuberance Induced Fluctuating Pressure Field, 
Upstream (0“ = 0°), 8-inch Diameter Protuberance, M = 1 .00, 

Re/ft = 3 x 10 6 00 
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Figure 53. Concluded 



o 

2.000 

o 

1.500 


r/u 

□ 

1 .875 

o 

1.375 

c 

1 .000 

A 

1.750 

A 

1.250 

□ 

0.750 

O 

1.625 

o 

1.125 

A 

0.625 


HI III 


i ■Him 


i ■hiui 


IUUI 


iBSnSSSSiSSSiSSSSiilSSSSiBSSiiSSSiiBB 


III 


IBM 


(■■I 


■■■El 


bt»H 

RhSw^;s5hiSS 
MMHBlftBaBMIMaHIE 

IIE!SB!SIII£l<IIIIBBII( 

IllIBB&SISliillIBBIII 

H B 

I iBBawtiMiil 

■■liimvii 

— — ■ | 

BMHiW WBBifwiil 

■MMflllB IB-lL-Jt\‘BBIEl 


mvm 


IRII 


■II 


■II 


jfcaefii 

■111 

II — ■■■«!■ 

■MMIIlilll 

1hb;i| 


IBIIIBBIIIBBIIIBBIII 

lllIBBHIIBBHIlfliilll 


BBIIWBIIIBBIIIBBIIE 

■bigsmiibbiiibEiii 

BBWBBBMWB W t E^JBB|toSlMryyi|y*^M jBEWIBB HIB ilBBBWBBMBH BC( 

BIlIBBIUSlIlIBBIIf 
BlllflfllllK’i^.IBBII! 

■III ■■ii!iil*5hhi 

■in 


■■■■i 


■re »B 5 il 

mm 


!BMB!!!BB!!liII!! 

IShhISBBSSiBBSBBSiSBmSBSe 

IBBBlIBBHilimiBIIBiiftlllE 


o. fi/°= 1.0 


Figure 54. Power Spectra of the Protuberance Induced Fluctuating Pressure Field 
Upstream (0" = 0°). 8-inch Diameter Protuberance, M = 1 .20, 
Re/ft = 3x10 6 °° 
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Figure 55. Power Spectra of the Protuberance Induced Fluctuating Pressure Field 
Upstream (0" = 0°), 8-inch Diameter Protuberance, M =1 .40, 
Re/ft = 3 x 10 6 00 





gs 

■■HHUHmHUUBBHHUUH 

bhihb mmtwmm ■■« 





Sun 

hhiihhb 

B 

■Si ISiiilBiiiii 

■IS 




■WIIKWIHHIHHII 


; i - .’**■■ ■■n ■■ wm m m m^m mmm w* n 

iSESiSSal^ ISSBBSi 


lX«£?ll 

j&aajs&M 


■■El 


■■■■ 


B — ■BIHMl'aAjBB M aHW 

immumiZimmmi 

ibiiibmbbi 


IHI 


ll'|III!ll 


iruKpiiwaiin 

IE&B&3IUUIII 

lil!B3&k!BBIIII 

■ wb «! w ^* i«i .var^i ■«■■■■ 

!B!SB£§ ii^B! 


Pill 

»,;>i 


m 


■■m 

Mill 


■■El 


llll 


miik^ 


■■El 


mi 


1.0 0.0001 


a. h/D = 1 .0 


Figure 56. Power Spectra of the Protuberance Induced Fluctuating Pressure Field 
Upstream (0~=O°), 8-inch Diameter Protuberance, M = 1.60, 
Re/ft = 3 x 10 6 °° 
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Figure 57. Power Spectra of the Protuberance Induced Fluctuating Pressure Field 
Downstream (0“ = 180°), 8-inch Diameter Protuberance, M = 0.60 
Re/ft = 3 x 10 6 00 
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Figure 60. Power Spectra of the Protuberance Induced Fluctuating Pressure Field 
Downstream (0" = 180°), 8-inch Diameter Protuberance, M =1 .20 
Re/ft = 3 x 10* °° 
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Figure 61 . Power Spectra of the Protuberance Induced Fluctuating Pressure Field 
Downstream (0“ = 180°), 8-inch Diameter Protuberance, M = 1 .40 
Re/ft = 3xl0 6 00 
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a . h/D = 1 .0, r/D = -1 .375 to -1 .250 


Figure 63 . 


Longitudinal Cross-Spectra for the Protuberance Induced Fluctuating Pressure Field 
Upstream (0 - = 0°), 8-inch Diameter Protuberance, M =1.20, Re/ft = 3x 10 6 , 
£/D = 0.125 00 
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Figure 64. Longitudinal Cross-Spectra for the Protuberance Induced Fluctuating Pressure Field, 
Upstream (9 = 0°), 8-inch Diameter Protuberance, M = 1 .40, Re/ft = 3 x 10 6 , 
I/D = 0.125 00 
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b. h/D = 2.0, r/D = -1 .875 to -1 .750 
Figure 64. Concluded 
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b . h/D = 2 .0, r/D = -1 .875 to -1 .750 


Figure 65. Concluded 
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Figure 66. Longitudinal Cross-Spectra for the Protuberance Induced Fluctuating Pressure Field, 
Downstream (0 - = 180°), 8-inch Diameter Protuberance, M =1.20, Re/ft = 3xl0 6 , 
£/D = 0.125 00 
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b . h/D = 2 .0, r/D = 1 .500 to 1 .625 
Figure 66. Concluded 
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a . h/D = 1 .0, r/D = 1 .500 to 1 .625 


Figure 67. Longitudinal Cross-Spectra for the Protuberance Induced Fluctuating Pressure Field, 
Downstream (0“- 180°) 8-inch Diameter Protuberance, M =1 .40,Re/ft=3 x 10 6 , 
£/D = 0.125 00 
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b. h/D = 2 .0, r/D = 1 .500 to 1 .625 
Figure 67. Concluded 
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Figure 68. Longitudinal Cross-Spectra for the Protuberance Induced Fluctuating Pressure Field, 
Downstream (0" = 180°), 8-inch Diameter Protuberance, M = 1 .60, Re/ft = 3 x 10 6 , 
£/D = 0.125 00 
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Figure 68. Concluded 
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a. h/D = 1.0, r/D= -1.375 to -1.250 

Figure 69. Longitudinal Narrow Band Convection Velocity for Protuberance Induced 
Fluctuating Pressure Field, Upstream (0"=O°), 3-Inch Diameter 
Protuberance, Mqq = 1 .20, Re/ft = 3xl0*, £/D = 0.125 
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a. b/D = 1.0, r/D - -1 .625 to -1 .500 

Figure 70. Longitudinal Narrow Band Convection Velocity for Protuberance Induced 
Fluctuating Pressure Field, Upstream (9" =0°), 8-Inch Diameter 
Protuberance, = 1 .40, Re/ft = 3 x 10 6 , I/D =0.125 
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a. h/D= 1.0, r/D = -1.375 to -1.250 


Figure 71 . Longitudinal Narrow Band Convection Velocity for Protuberance Induced 
Fluctuating Pressure Field, Upstream (0" =0°), 8-Inch Diameter 
Protuberance, Mqq = 1 .60, Re/ft = 3 x 10 6 , k/D = 0.125 
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b. h/D = 2. 0, r/D = -1.875 to -1.750 
Figure 71 . Concluded 
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a . h/D = 1 .0, r/D = 1 .500 to 1 .625 

Figure 72. Longitudinal Narrow Band Convection Velocity for Protuberance Induced 
Fluctuating Pressure Field, Downstream (0 = 180°), 8-Inch Diameter 
Protuberance, - 1 .20, Re/ft = 3xl0 6 , l/D = 0.125 
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a. h/D = 1 .0, r/D = 1 ,500 to 1 .625 

Figure 73. Longitudinal Narrow Band Convection Velocity for Protuberance Induced 
Fluctuating Pressure Field, Downstream (0“ = 180°), 8-Inch Diameter 
Protuberance, = 1 .40, Re/ ft = 3 x 10 6 , £/D = 0.125 
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a. h/ D = 1 .0, r/D = 1 .500 to 1 .625 

Figure 74. Longitudinal Narrow Band Convection Velocity for Protuberance Induced 
Fluctuating Pressure Field, Downstream (0’ = 180°), 8-Inch Diameter 
Protuberance, Re/ft = 3 x 10 6 , l/D = 0.125 
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Figure 75 . Composite Schematic of Protuberance Flow Field Characteristics, 
M = 1 „60, h/D = 2.0 
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a. Side View 


Figu.e 76. Titanium Dioxide Flow Patterns Over the Wall of Cylindrical Protuberance 

















Figure 77. One-Third Octave Spectra for Typical Microphone Locations Upstream of the 
Protuberance, Zero Level = 120 dB, M = 1 .60, h/D - 2.0 
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b. r/D = -2.125 to -2.000, 4/D = 0.125 



Figure 78. Variation of Longitudinal Cross-Spectra Within the Upstream Separated 

Flow Field, M = 1 .60, h/O = 2.0 
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Figure 78. Concluded 
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Figure 79. 


Smoke Filament Flow Pattern Upstream of a Cylindrical Protuberance 
at Low Speed (Reference 39) 
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